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National Bureau of Standards 
Materials Measurements 

Summary 

This report describes NBS work for NASA in support of NASA's Materials 
Processing in Space Program under NASA Government Order H-27954B (Properties 
of Electronic Materials) covering the period April 2, 1979 to April 1, 1980. 

The work emphasizes .aterials measurements an’ has two main thrusts: 

1) Carrying - at precision measurements in space and investigating the 
feasibility of improved measurements when the space environment offers a 
unique opportunity for performing such measurements. These measurements 
would be useful for either space processing or processes on the ground. 

2) Obtaining precision measurements on materials properties when these 
properties are important to the design and interpretation of space processing 
experiments. These measurements would be carried out either in space or on 
the ground. 

This work has been carried out in four tasks. These tasks have, as two 
of their focal points, the role of convection effects and the role of container 
effects, both of which would differ in space-based experiments from those 
found in ground-based experiments. The results obtained for each task are 
given in detailed summaries in the body of the report. Briefly, in Task 1 - 
Surface Tensions and Their Variations with Temperature and Impurities - 
measurements in an Auger spectrometer of surface impurity concentrations on 
liquid gallium showed that the principal impurities were oxygen ard carbon. 

The impurities showed a tendency to collect into plates or clumps. Heating 
in hydrogen atmospheres and Ion sputtering were both found to be effective 
methods of removing impurities. When these methods of removing impurities 
were applied, reproducible surface tension values at room temperature 
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were found with y = 710 mJ/m . The surface tension decreased quadrat ically 

2 

with increasing temperature to y = 690 mJ/m at 500 °C. In Task 2 - Solutal 

Convection During Directional Solidification - in Pb-rich Pb-Sn off-eutectic 

alloys, macrosegregation caused by solutal convection was not reduced by 

vertical or horizontal field*: of 0.1T, but downward solidification (liquid 

below sclid) virtually eliminated macrosegregation in small 3 mm) diameter 

samples. Calculations were made to predict the onset of convective and 

interfacial instabilities during directional solidification of lead-tin 

alloys including effects of variations in gravitational acceleration and 

Prandtl number. Results obtained in calculations with a constant gravitational 

acceleration of 10 ^ g g show a noticeably larger stability range than those 
-4 -6 

for 10 g g ; at 10 g g the instabilities will develop more slowly with tine. 

In Task 3 - A Thermochemical Study of Corrosive Reactions in Oxide Materials - 

phase assemblages of selected compositions on the joins RfFe^ .-SIq 5 ] 0 ^ - 

SiO^ and KFeO^ - were determined over a large range of oxygen partial 

pressures and the temperature range 800 °C to 1400 °C. In melts containing 

ferrous ions, excessive reaction with platinum containers occurred. In Task 

4 - Thermodynamic Properties of Refractory Materials at High Temperatures - 

use of pyroelectric detectors to determine the radiant heat loss from 

spherical samples as cooling occurs in free-cooling experiments is being 

investigated. The heat capacities and, under some conditions, the thermal 

conductivities in leviatated high temperature materials can be derived from 

simultaneous measurement of radiant heat loss and rate of temperature change. 

Details of the thermal gradients to be expected in both metallic and non- 

metallic samples have been examined. Limitations from sample vaporization 

-4 

will occur when the vapor pressure reaches the vicinity of 10 Pa (10 atm). 
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Task 1 


Surface Tensions and Their Variations 
with Temperature and Impurities 

S. C. Hardy 
Metallurgy Division 
Center for Materials Science 

and 

J. Fine 

Surface Science Division 
Center for Thermodynamics and Molecular Science 

Summary 

The surface tension of gallium has been measured as a function of tem- 
perature in vacuum and hydrogen atmospheres using cups of quartz and graphite 
to contain the sessile drops. The results are generally in good agreement 
with two previous measurements and show that the surface tension decrease is 
nearly quadratic in temperature. A non-linear least squares analysis gave 
the following function as the best fit to the data; 

Y - 709.9 + .0016 (T-29.87) -.0000923 (T-29.37) 2 . 

Measurements in vacuum using quartz cup-, however, have revealed a complex 
heating history dependence which we think is due to an oxide contamination of 
the gallium by the quartz. 

In related work the surfaces of liquid gallium sessile drops have been 
studied ir an Auger spectrometer to determine the concentration of impurities. 
The gallium surface after formation of the drop was found to be covered with 
a solid layer containing a high concentration of oxygen and carbon. Prolonged 
heating to several hundred degrees effectively desorbed the oxygen but the 
carbon was still present in the form of a precipitate. The conditions under 
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which these precipitates form is not understood. The surfaces were also 
cleaned by sputtering with argon ions and complex flow phenomena observed. 
The surface tension at the melting point of a drop with a sputtered clean 
surface was in fair agreement with the results of the sessile drop measure- 
ments. 
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Introduction 


Studies of liquid gallium surfaces have proceeded along two lines in 
the past year. Extensive surface tension measurements have been made using 
the apparatus and techniques described in last year's report. In addition, 
we have . tempted to characterize the surface composition of liquid gallium 
in an Auger spectrometer. These latter experiments have revealed complex 
impurity phenomena which we do not fully understand. 

The central prooxem in -^cermining the surface tension of liquids is 
to produce a clean surface. Segregation of impurities at the liquid-vapor 
interface will drastically plter the surface tension as described by the 
Gibb's equation. This sensitivity to surface active impurities is probably 
the basic cause for the wide discrepancies in reported surface tension 
values. A related problem is that there is almost no independent character- 
ization of the chemical composition at the interface with a surface tension 
measurement. It is assumed that the procedures and materials used have 
produced a clear surface, the cleanest surface being ascribed to the one 
with the highest surface tension, i.e. , the one with the least surface 
segregation. An independent characterization of the surface by some 
direct technique would be extremely valuable in the interpretation of 
surface tension measurements, particularly in cases where there is a wide 
discrepancy in the results. The Auger work was undertaken primarily to 
address this problem. 

Oxygen is probably the most troublesome surface contaminant for most 
metals; gallium, whi< is chemically similar to alumirmm, is particularly 
subject to the formation of stable surface oxides. These may be either 


5 


Ga^O or aGa^O^^^. There are four additional modifications of the trioxide 
which are formed only under unusual conditions. Several approaches are 
available which one hopes might succeed in producing a gallium surface free 
of oxide. Our samples are prepared in a syringe under alcohol containing 
HC1 which dissolves gallium oxides. The gallium is squeezed out into a 
shallow cup to form the sessile drop under an atmosphere of flowing nitrogen 
or helium. Although gross oxides are removed by the acid, a chin oxide 
film undoubtedly remains after this procedure. Observations supporting 
this will be discussed later in connection with the Auger work. Once the 
sessile drop is formed in the vacuum system, t may be cleaned farther Hv 
heating. This desorbs weakly bound gases an . at around 500 °C, dec es 
the trioxide to the suboxide according to the following reaction. 

Ga^O-j + 4Ga +~ 3Ga 2 0 

The suboxide has a high vapor pressure and is deposited on the walls of tne 
vacuum system. This reaction, however, does not go to completion as shown 
by x-ray diffraction measurements^. Alternatively, the surface oxide can 
be reduced by heating in hydrogen by the following reactions: 

Ga 2 0 3 + 2H 2 t Ga 2 0 + 2H 2 0 

Ga 2 0 + H 2 * 2Ga + H 2 0 

We have used both of these approaches in attempts to produce clean surfaces. 
A third possibility for cleaning the surface is ion '.c sputtering. We will 
discuss this technique in connection with our Auger work. This alternative 
is very attractive because it can be carried out at low temperatures, thus 
minimizing the diffusion of gases and other impurities from the cup into 
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the gallium and any reaction of the gallium w'.th the cup or residual gases. 

A detailed discussion of "he apparatus and techniques used in this 

[ 2 ] 

work is contained in last year's report . Rather chan repeat this, we 
will simply note that the gallium drops are contained in cups of various 
materials in an ion pumped vacuum system capable of reaching the 10 ^ 
torr range. The optics and film used in the measurements have a resolution 
which permits the surface tension to be measured over the profile of the 
drop with a precision of about 0.5% using the Bashforth and Adam 1 technique. 
This was verified last year with a computer program which fitted the entire 
measured profile tc calculated solutions of the Young- cap lace equation. 

Thus we have some confidence in the measurements themselves. Additional 
experimental details will be given when necessary with the presentation of 
the data. 

Surface Tension Measurements 

The initial surface tension measurements used quartz cups to contain 

the sessile drops. Fig. 1 reproduces data shown in last year's report. 

Although the surface tension appears to decrease linearly with temperature 

above 120 °C, there is a maximum in the curve at that temperature. This 

kind of behavior is usually interpreted as evidence of impurity segregation 

which builds up at low temperatures because desorption is too slow to keep 

the surface clean. This seemed to be a reasonable explanation because 

hours were required to cool the drop to this temperature range after the 

initial high temperature bake, a time sufficient to adsorb many monolayers 

of gas from the residual atmosphere of the vacuum system which was at a 
_ & 

pressure of 10 torr. 
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Analysis of subsequent measurements proved this adsorption explanation 

to be wrong in that impurity did not originate in the residual gases in 

the vacuum system. When the drop was heated for twenty hours at 570 °C and 

then cooled, the data shown as open circles were obtained. The lengthy heat 

treatment enhanced the low temperature depression of the surface tension and 

moved the maximum in the data to a higher temperature. An additional twenty 

hour heat treatment lowered the surface tension ever, more dramatically as 

shown by the open boxes. The vacuum level during these measurements was 
—8 

about 10 torr. In subsequent experiments on other drops the entire vacuum 
system was baked and the pressure was in the mid 10 ^ torr range. Measure- 
ments at this lower pressure confirmed the behavior shown in Fig. 1, 
suggesting the residual gases were not the adsorbed spec’es. 

This observation that the temperature and heating time control the 
behavior at lc..er .mperatures indicates a contamination of the sample by the 

crucible at high temperatures. Althougn gallium is known to attack quartz at 

[3] 

temperatures over 1000 °C , it was thought to be non-reactive at ^he 
maximum temperatures used in our experiments. It is possible, however, that 
the contamination does not come about from gallium dissolving the crucible 
but rather from gases or other components of the crucible diffusing into the 
gallium. Alternatively, t^e silica glass tube forming the vacuum jacket 
could release gases during bakeout which react at the gallium surface. A 
likely candidate for this would be water which is the primary constitutent 
desorbed fron quartz at 600 °C L J . Thus a possible reaction taking place 
might be the following: 
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2Ga + H 2 0 * Ga 2 0 + H 2 

Another possibility for contamination exists, however. Silicon contamination 

of gallium in quartz containers has been tudied in connection with the growth 
[41 

of GaAs 1 

Support for the oxide contamination hypothesis is shown in Fig. 2. This 
sample was loaded into the cup as a solid under high vacuum conditions in an 
attempt to avoid adsorbed gases at the quartz-gallium interface. The gallium 
was frozen under alcohol in the form of a cylindrical ingot with a thin 
rhenium wire along the axis. The rhenium ware was attached to iron slugs 
which could be moved magnetically through the silica wall of the vacuum 
system. With this solid loading technique, the cup could be baked at high 
temperature and the gallium loaded into it without breaking vacuum and back 
filling with an inert gas, a procedure which undoubtedly results in some 
oxygen and water adsorption on the quartz cup. When this procedure v.'as 
followed and the gallium subsequently heated to 760 °C for twenty-four hours, 
the surface tension varied with temperature as shown by the solid points in 
figure 2; a very strong surface tension depression at low temperatures with a 
maximum at about 300 °C. This behavior resembles that shown in Fig. 1. 

Hydrogen gas was then admitted to the tub^ through a liquid nitrogen cooled 
trap to about atmospheric pressure. When the gallium was heated to 550 °C for 
30 minutes, the surface tension at low temperatures increased dramatically and 
effectively removed the maximum in the curve. These data are shown as open 
circles in Fig. 2 and suggest that an oxide at the surface is being reduced by 
heating in hydrogen, the second cleaning procedure mentioned in the Introduction. 
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The surface oxides in this experiment may not have been Identical to 

those formed with the drop of Fig. 1 because of the different procedures 

used to form the sessile drop. There was probably a heavy oxide layer 

present when the gallium was loaded as a solid because the ingot was 

exposed to air during transfer to the vacuum system. Larger surface 

tension depressions were found at low temperature in this case than when 

the liquid loading procedure was used. This high concentration of oxide 

may have led to silicon contamination of the gallium by dissolving the 

quartz during the extensive high temperature bake. 

A third experiment was performed in the quartz cup. Prior to forming 

the sessile drop the system was evacuated and the cup and adjacent tube 

baked at 900 °C for eighteen hours in an attempt to thoroughly eliminate 

water. The gallium was then added to the cup as a liquid and the system 

_8 

evacuated to the 10 tore range. After heating the gallium to about 

100 °C to desorb weakly bound gases, hydrogen was admitted to the tube 

through a copper coil immersed in liquid N£- Heating to progressively 

higher temperatures for periods o f about 30 minutes and then cooling to the 

melting point and measuring the surface tension established that significant 

cleaning of the surface occurred at temperatures as low as 316 °C. After 

30 minutes at 460 °C, the low temperature depression of the surface tension 

was found below 50 °C and was so slight it could not positively be separated 

2 

from measurement imprecision, i.e. , it was about 7 mJ/m , which is 1Z of 
the surface tension. The gallium was then heated to 550 °C for fifteen 
minutes and photographed. The surface tension was measured after equilibration 
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at a number of lower temperatures. These data are shown as filled circles 
in Fig. 3. We also show as filled squares the results obtained after 
heating, to 600 °C for eighteen hours. This hotter and longer heat treatment 
seems to have raised the surface tension at the melting point, but the 
change is only 2 mj/m^, well within the normal measurement imprecision. The 
general agreement of the two sets of data is good. 

We see, therefore, that prolonged heating in hydrogen with the 

quartz cup seems to have little effect on the surface tension at lower 

temperatures in contrast to the behavior in the vacuum experiments typified 

in Fig. 1. Moreover, the surface tension now varies monotonically with 

temperature; there is no identifiable relative maximum in th' data. The 

temperature variation, however, does not seem to be linear. The departure 

from linearity could be caused by impurities in the hydrogen gas being 

adsorbed at low temperatures. An increase in surface tension of 1% to about 
2 

718 mJ/m at the melting point would effectively linearize the data. 

After the measurements described above, the system was evacuated. When 

-8 

the pressure had dropped to the high 10 torr range, the gallium was heated 
for about 15 hours at 550 °C. A subsequent measurement of surface tension 
as a function of temperature produced the data plotted as open triangles in 
Fife. 3. The surface tension has been reduced significantly, but the 
variation with temperature is still monotonic. This is quite different from 
the results cited earlier and illustrated in Fig. 1. It is not clear 
whether the change has been brought about by the extended high temperature 
treatment prior to forming the drop or by initially heating in hydrogen 
rather than vacuum. The results using quartz cups seem to be very 
dependent on the history of the sample; i.e., the atmosphere and pr ivious 
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heat treatment. Although the results in hydrogen are encouraging, the 
depression of surface tension with heating in vacuum indicates a container 
reaction which raises questions about the sample purity in all of the 
experiments. Presumably the contaminant is a 11 oxide which is reduced by 
the hydrogen, but there may be an associated silicon contaminant also 
playing a role. 

The ambiguous results with quartz led us to try a number of otner cup 
materials. Sessile drop experiments were performed in cups of alumina, 
boron nitride, and two types of graphite, for the most part these materials 
were inferior to the quartz. Cups of high purity, pyrolitic graphite, 
however, gave more reproducible and consistent results than the other 
materials. A series of measurements were made under a variety of conditions 
; one of these graphite* cups in an attempt to see if the surface tension 
changed systematically with heat treatment and atmosphere. To eliminate 
considerations of changing drop sizes and cup orientations or uncontrollable 
differences in the gallium introduced during preparation, all of the data 
we 'ill show have been taken with the same drop. These results, however, 
are reproducible from drop to drop. 

After the drop was formed using the liquid loading technique, the 

system pressure quickly reached the low 10 ^ torr region. The drop was 

heat id overnight to about 85 °C to desorb weakly attached gases. The 

2 

surface tension was under 600 mJ/m at this temperature afterwards in- 
dicating no significant cleaning of the surface had occurred. The temperature 
was then raised to 185 °C for two hours, the pressure in the system being 
about 10 ^ torr. When the drop was cooled to the melting point, the surface 
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tension was measured to be 696 mJ/m indicating a substantial cleaning of 

the surface had been achieved. Subsequent heating for short periods to 

progressively higher temperatures slightly increased the melting point 

surface tension. Fig. 4 shows that the temperature variation of the 

surface tension after heating to 465 °C is approximately linear with no 

evidence of a maximum. The vacuum level during these measurements was 
—8 

about 2 x 10 torr. Subsequently the drop was heated to 565 °C for 
eighteen hours with the results shown in Fig. 5. These results do not 
differ significantly from those shown in Fig. 4. The Important fact is 
that the measured surface tension has not been changed with a prolonged 
high temperature bake, in contrast to the results with quartz cups in 
similar vacuums. 

After these measurements hydrogen was admitted to the vacuum system in 
the normal manner and the temperature was raised briefly to 550 °C. This 
produced ro significant changes from the data of Fig. 5. An overnight 
'.aeating to 600 "C produced the results shown in Fig. 6. Again, there does 
not seem to be any significant change from the results in vacuum shown in 
Fig. 5. 

A direct comparison of the data obtained after prolonged high temperature 
heating with graphite cups in vacuum and hydrogen and with quartz cups in 
hydrogen is made in Fig. 7. We have simply plotted the data shown in Figs. 

3, 5, and 6 on one graph; no adjustment or normalization has been used. As 
is evident, the measurements are in agreement and suggest that the surface 
tension of gallium does not vary linearly with temperature. A least squares 
non-linear fit to an equation of the form 
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Y - a - b (T-29.87) - c (T-29.87) 2 
gave final parameter estimates of 

a - 709.92 (.57) 

b - -.160275 x 10 -2 (.5311 x lO* 2 ) 
c - .923122 x 10' 4 (.9474 x 10~ 5 ) 

where the numbers in parentheses are the approximate standard deviations. For 
b the standard deviation is larger than the estimate of the parameter indicating 
that this coefficient is negligible. Thus, the dependence is essentially 
quadrat ic . 

The general agreement of the vacuum and hydrogen measurements in graphite 

with the hydrogen measurements in quartz reinforces our feeling that gallium 

picks up an oxide impurity when it is heated to high temperature in quartz. 

This oxide presumably is effectively reduced in hydrogen but in vacuum it 

segregates at the surface and depresses the surface tension at low temperatures. 

A superposition of Fig. 1 on Fig. 7, however, reveals an additional complication. 

Although the high temperature data of Figs. 1 and 7 superpose nicely, at and 

below 200 °C, the surface tension in the quartz cup in vacuum rises above the 

data of Fig. 7. The maximum of Fig. 1 which occurs at about 150 °C is about 6 
2 

mj/m higher than the curve of Fig. 7. This is not a large discrepancy but a 
reexamination of the original data Indicates that it is real. It may result 
from the fact that gallium oxides seem to form in patches rather than in a 
uniform surface layer. Thus, it is possible that the surface is no longer 
strictly described by the Young-Laplace equation and the Bashforth and Adams 
analysis is invalid. 
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Auger Spectrometer Experiments 


A 8 mentioned previously, the lack of any independent characterization of 
the surface is a serious problem in surface tension measurements. A number of 
modern surface analysis techniques are available which. might fill this need, 
but they have been applied infrequently to liquids. The high vapor pressure 
of most liquids, of course, reduces the range of materials which can be studied. 
Nevertheless, the complementary use of these techniques with surface tension 
measurements seems feasible for some liquids and such studies are being 
attempted by several groups. 

Liquid gallium is compatible with high vacuum instruments because of its 
low vapor pressure and melting point. It melts near 30 °C and the vapor 
pressure is about 10 ^ torr at 700 °C so that a large temperature range is 
usable. We have been studying liquid and solid gallium in an Auger spectrometer. 
Our specimens are in the form of sessile drops and are contained in cups of 
graphite or rhenium which can be heated to over 500 °C and cooled to low 
temperatures. The drops are formed using the liquid loading technique described 
previously. The geometry of the spectrometer is not ideal for these measurements 
because the drops must be positioned well off the axis of the hell jar in 
order to be within working range of the horizontal cylindrical mirror analyzer. 
This makes simultaneous Auger analysis of the surface and photography of the 
drop profile impossible. Furthermore, the cup cannot be precisely aligned 
with the horizontal. Thus surface tension measurements in the present spectro- 
meter are less accurate than those we have been discussing. 

Although the surface of the drop during Auger analysis is difficult to 
observe visually, a video image can be obtained by a technique similar to that 
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used In scanning electron microscopy (SEM). A primary electron beam Is 
raster ed over the surface of the specimen and the local change in secondary 
electron emission is used to Intensity modulate a video monitor in synchronization 
with the two-dimensional position of the electron beam. This Imaging mechanism 
reveals surface features which are difficult to see optically. Although the 
video Images do not have high contrast or resolution, they provide valuable 
information about the sunace. Furthermore, the video imaging permits the 
sample orientation to be adlusted so that selected regions may be analysed 
with the Auger spectrometer. The instrument is equipped with an ion gun and 
a gas handling system so that the sample surface can be sputtered. In 
addition, a quad ru pole mass spectrometer permits residual gas analysis. The 
pumping is provided by liquid nitrogen cooled sorption pumps and ion pumps 
so that the system is oil free. 

After forming the sessile drop in the cup and evacuating the system to 
_9 

the low 10 torr region, an Auger spectrum of the surface has intense 

oxygen and carbon lines as seen in Fig. 8. Although the video Image of the 

drop is uniform and featureless as would be expected of a liquid, the 

surface is actually solid. If the electron gun is left on for a number of 

hours, a dark rectangular spot develops where the beam has been rastered 
over the surface which persists unchanged and stationary after the beam is 
turned off. The amplitude of the Auger lines vary differently over this 
dark region, oxygen decreasing and carbon Increasing as the beam is moved 
from the center of the raster area to an unbombarded area. This does not 
necessarily mean that the concentration of these constitutents is varying in 
the same way because the Auger line intensities of different elements are 
interrelated. 
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The loading technique, therefore, produces a sessile drop with a solid 
oxide surface layer containing significant concentrations of carbon. Our 
experience with the surface tension measurements suggested this layer might 
be decomposed and desorbed by heating to a moderate temperature. The effect 
of heat treatment on the gallium surface was examined by heating the drop 
briefly to progressively higher temperatures, then cooling and taking Auger 
spectra. Aftar heating to 250 °C for five minutes, the oxygen line was 
halved and the carbon line was also reduced somewhat. However, the surface 
was no longer uniform and featureless. A "white" streak (as seen on the 
video picture) developed which gave Auger lines for oxygen and carbon at 
least twice the amplitude of the adjacent, featureless surface. Thus it 
appears that the thermal decomposition and desorption of the surface layer 
is not uniform, or that segregation and growth may occur during the heating 
and cooling procedure. The surface of the drop after heating is, of course, 
no longer solid and flows can develop. As the drop cools a phenomena is 
seen which we do not yet fully understand. The drop surface is initially 
featureless after heating to 400 °C. On cooling to about 120 °C, however, 
small black spots appear near the apex of the drop and move down the side 
collecting into a broad band on the surface above the cup edge. An Auger 
spectrum of the deposit shows carbon but no noxygen. This process was 
observed in both rhenium and graphite cups. The Auger spectra of areas of 
the surface without precipitates have oxygen and carbon lines which have 
diminished to the point at which they are indistinguishable from the noise. 
Concurrently the gallium lines become more intense. 
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Other experiments have shown tteit heating to 200 °C for several hours 

produces a surface which Is clean except for the precipitates described 

above. The surface of the precipitates remains clean for long periods; no 

oxygen or carbon was detected after twenty-four hours at a pressure of 6 x 
—8 

10 torr. Exposure to the residual atmosphere for a number of days produced 
a significant carbon line and a somewhat smaller oxygen line. This is to be 
expected because CO is an important residual gas in the vacuum system 
presumably because of the hot fU ,-ents in the ionization gauge. 

In addition to carbon and oxygen the measurements in the vitreous 
graphite cup found some precipitates with a high sulfur content. An Auger 
analysis of the cup wall also found some sulfur. We suspect that the 
sulfur was introduced in a grinding procedure used on the cup during modi- 
fication. A more stringent cup cleaning procedure eliminated or greatly 
reduced sulfur indications in later experiments. Several precipitates on 
one drop gave strong lead lines. 

Another surface cleaning procedure we have tried is ionic sputtering. 
Within ten seconds of the initiation of sputtering with 2.5 keV argon ions, 
the uniform, solid surface film breaks up into large islands which move and 
rotate about the bombardment area. These islands shrink in the ion beam and 
eventually, when the surface is sufficiently clean, fluid flow develops 
which brings small islands from the rest of the drop to the site of sputtering. 
When the ion beam is turned off, the motion of the islands abruptly halts. 

With continued sputtering, the number and size of the islands decxeases. 
Eventually the flow slows and stops with a few residual islands resting near 
the cup ^dge below the ion beam impact area. An Auger spectrum of the surface 
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after sputtering for several hours is shown In Fig. 9. This Is identical 
to the result found after thermal cleaning and shows no carbon or oxygen. 

The same spectrum is obtained if the gallium is frozen and sputtered further. 

Although there are other explanations for the processes seen on the 
surface during sputtering, the generation of fluid flow by surface tension 
gradients comes immediately to mind. In this model the sputtering produces 
and maintains a clean area of high surface tension. If adjacent areas have 
a lower surface tension because of adsorption, a high surface tension 
gradient will exist about the periphery of the sputtering site. The 
tangential shear stress generated on the liquid surface by the gradient 
will set up a fluid flow on the surface which will be directed from the 
area of low surface tension into the area of high surface tension, i.e., 
into the sputtering site. Presumably this becomes a convection node with 
the surface flow entering the drop volume at that point. If the sputtering 
is interrupted, the surface quickly equilibrate' and the flow stops. As 
the surface becomes cleaner, the gradient is reduced and the fluid flow 
disappears. 

A surface cleaned by sputtering will apparently stay clean for days 
with no change in the Auger spectrum or any additional island formation. 

If, however, the gallium is heated briefly, on cooling a precipitation 
process occurs which is similar to the one described previously. Small 
dark islands appear and descend from the apex of the drop piling up in a 
belt of precipitates at the base. An Auger spectrum of these precipitates 
again shows a strong carbon line and no oxygen. Host of these precipitates 
disappear under ion sputtering within a minute which indicates they are 
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very thin. Repeated sputtering and heating seems to reduce their numbers, 
but it is not clear that they can be totally eliminated in this way. 

Our understanding of these sputtering and precipitation phenomena is 
limited with major questions to resolve about the origin and evolution of 
carbon in the gallium. It does not seem likely tha carbon concentration 

can be only that originally adsorbed on the surface, ably as C0 0 . We 
suspect a continuing source of contamination which is logically the graphite 
cup. Occasionally large and thick carbon islands have been seen which we 
feel must liave originated in the cup. Thus it is possible that the cup 
surface is flaking off under repeated heating and cooling cycles. It does 
seem that there was significantly l^ss carbon with rhenium cups. However, 
this work is very limited because the gallium drops tend to run out of 
rhenium cups probably because of good wetting. Some carbon is certainly 
due to long exposures to the residual gases in the vacuum system because it 
is uniformly present on the surface and not associated with a precipitate, 
at least before heating. It is possible that the sputtering does not 
effectively remove carbon from the system but just disperses it into the 
vu..jme as very fine particles so that it is not seen. Heating then could 
lead to agglomeration and reappearance at the surface because of buoyancy 
and fluid flow. 

As mentioned earlier we cannot make accurate surface tension measure- 
ments on the sessile drops in the Auger spectrometer because of geometric 
and alignment, problems. We have, however, measured the surface tension of 
a drop which was sputtered for several liours and whose Auger spectrum 
indicated « cl^an surface. In order to make the measurement, the drop was 
repositioned in the Auger spectrometer after sputtering and analysis so 
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that its profile could be photographed in parallel light. The surface 

2 

tension was about 700 mJ/m . Recogni;. Lng that this may t i inaccurate, it 

seems unlikely that the error could be more than a few percent. It seems 

2 

beyond possibility that the surface tension could be as high as 800 mJ/m , 

[ 2 ] 

the anomously high values found previously in two instances . A new 
spectrometer is being assembled which will permit accurate surface tension 
measurements and simultaneous Auger spectral analysis. We believe that it 
will permit a resolution to this discrepancy in previous measurements. 

Discussion 

The experiments with gallium in the Auger spectrometer have shown the 
primary surface impurities to be oxygen and carbon. A relatively low 
temperature heating in vacuum can clean the surface of oxygen but the carbon 
is either not effectively desorbed or is being resupplied, possibly from the 
crucible. The carbon after heating is contained in thin plates which 
collect on the drop surface near the cup edge. Thus the surface of the 
drop, although clean for the most part, will have regions in which it is not 
liquid. The entire drop shape in this case is not described precisely by 
the Young-Laplace equation. Because the precipitates collect at the base, 
however, the resultant distortion of the drop near the apex may be small. 
Thus the surface tension values calculated by the Bashforth and Adams 
technique above the d >p maximum may not be seriously affected. Assuming 
such a precipitation phenomenon is also occurring in our sessile drop 
experiments, the fact thac very nearly the same surface tension is cal- 
culated at various points on the drop in this region suggests the profile 
is not severely perturbed. We do not know, of course, that precipitation 
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is occurring in the sessile drop experiments. The cups are different forms 
of graphite from different sources in the two experiments; the Auger cup is 
"vitreous" graphite and the sessile drop cup is high purity pyrolitic 
graphite. Sessile drop measurements in a companion vitreous graphite cup 
indicated some surface tension depression at low temperatures, but it is 
dif icult to compare the results because of variations in procedure in the 
two groups of measu aments. The approximate agreement at the melting point 
of the surface tension values shown in Fig. 7 with that measured for a drop 
that was sputtered clean in the Auger spectrometer is additional evidence 
that the precipitates, if present, are not having much effect on the drop 
shape. 

The sessile drop surface tension measurements and the measurements and 
observations in the Auger spectrometer demonstrate a strong sensitivity on 
cup material and heating nistory. Certainly quartz cups used at high 
temperatures in vacuum result in contamination of the gallium with what 
seems to be an oxide because of the ease with which it is reduced in 
hydrogen. The gallium may also be contaminated wich silicon which has 
little effect on surface tension values presumably because of low con- 
centration or no tendency for surface segregation. 

The surface contamination by the residual gases in the vacuum system 

is much less significant than that originating in the cup. We have found 

in the sessile drop experiments that surface tension values at the melting 

—8 —7 

point decline very s^wly i vacuum levels of 10 to 10 torr. This is 
also true in the hydrogen atmospheres. The Auger measurements support, this 
observation. 
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The general agreement of the surface tension values measured in 
graphite cups in vacuum and hydrogen with the values found in quartz cups in 
hydrogen is shown in Fig. 7. Although this is strong evidence the gallium 
surfaces are in a similar state for these three groups of measurements, it 
does not prove the surfaces are clean. The very limited measurements of 
surface tension on the one drop in the Auger spectrometer, however, gave 
surface tensions at the melting point which are only about 1Z lower than 
those shown in Fig. 7. Given the probable error in these measurements, 
this is good agreement and suggests that the data of Fig. 7 are typical of a 
clean surface. 

There are, however, unanswered questions and unresolved contradictions 

in the surface tension measurements. Fig. 1 is data taken in a quartz cup 

which shows a low temperature surface tension depression which we have 

interpreted as arising from surface segregation of an oxide contaminant. A 

comparison of this data with the compilation of data from "clean" surfaces 

in Fig. 7 shows a very good agreement at higher temperatures. However, 

between 250 °C and the melting point, the surface tensions in the oxide 

contaminated gallium of Fig. 1 rise above that for the clean gallium by as 
2 

much as 6 mJ/m at the maximum. Ue, therefore, have a contradiction: if 

surface segregation is producing the maximum in Fig. 1 by depressing the 
surface tension at low temperatures, how can the surface tension of the 
contaminated surface be higher than that of a clean surface? A possible 
resolution would be that an impurity is also active in the data of Fig. 7 
and the curvature evident at low temperatures is due to adsorption. This 
impurity would necessarily be independent of cup material and atmosphere. 
Another hypothesis is that the impurities in the "contaminated" gallium 
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are confined to patches so that the drop shape is not being described by 
the Young-Laplace equation. We have no evidence supporting either of 
these speculations. 

As mentioned in last year's report, there have been three previous 
measurements of the surface tension of gallium as a function of temperature. 
The functions suggested by the authors to describe their data are as 
follows: 

Yj - 708 - .0031(T-29.87) - .00006>(T-29.87) 2t5 * 

Y XI = 706.7 - . 000647 (T-29. 87) - . 0000965(1-29. 87) 2 
Yin = 718.0 - .101(T-29.87) t?1 . 

The non-linear least squares analysis of the data shown in Fig. 7 mentioned 

previously gave the relationship 

Y = 709.9 + .0016 (T-29. 87) - . 0000923 (T-29. 87) 2 . 

Note that our results and those described by the first two relations above 

2 

differ at the melting point by no more than 3.2 mJ/m . Since small systematic 
errors must exist in all work, it is interesting to compare the various 
measurements by making them all coincide at one temperature. We do this 
in Fig. 10 where we have translated Yj and Yjj to agree with Fig. 7 at the 
melting point. Our data is then seen to lie between the two curves over 
the entire temperature range. The differences in the data adjusted in 
this way is hardly more than the precision with which the measurements 
were made in our work except possibly at the high temperature end of the 
curve. 

We also show the linear relationship (unadjusted). These 

measurements clearly disagree with the other data shown in Fig. 11. The 
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linearity of the surface tension values as a function of temperature in 

this work was cited as evidence that the gallium was cleaner than thnt. 

used in the measurements showning non-linear behavior. However, at 

temperatures above 125 °C, the surface tension values in the linear oata 

are lower than the non-linear data. At 360 °C, the highest temperature 

2 

used, they find a surface tension value about 15 m.T/m lower than the 
other results in Fig. 11. Thus it is hard to argue that the linear data 
represents a cleaner surface. This is the same problem discussed previously 
with respect to our own c'ata and we cannot rationalize the contradiction. 

We note, however, that the gallium preparation used in tnis work involved 
heating to 1000 °C in a quartz ampoule. We suspect their gallium may have 
been contaminated with silicon. We do not know how this would have affected 
their surface tension values. 
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Fig. 9. The Auger electron spectrum ol' a liquid gallium surface which was 
sputtered with 2.5 keV Argon ions for several hours. 
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The effect cf solute gradient induced convection during vertical 
solidification on the macrosegregation of Pb-rich Pb-Sn off-eutectic 
alloys is determined experimentally as a function of composition and 
growth rate. In many cases macrosegregation is sufficient to prevent the 
plane front solidification of the alloy. The transition from dendritic to 
composite structure is found to occur when the composition of the solid is 
close enough to the eutectic composition to satisfy a stability criterion 
based on G^/V (liquid temperature gradient/growth rate). A vertical or 
horizontal magnetic field of 0.1T (1 kilogauss) does not reduce macro- 
segregation, but downward solidification (liquid below solid) virtually 
eliminates macrosegregation in small {^3 nm) diameter samples. 

Additional calculations of the onset of convective and interfacial 
instabilities during the directional solidification of lead-tin alloys have 
been carried out with particular emphasis on constant gravitational 
accelerations of 980 (10 °) cm/s . 
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A. SOLUTAL CONVECTION INDUCED MACROSEGRECATIO:: AND Hit DENDRITIC TO 
COMPOSITE TRANS ITION IN OFF-EUTECTIC ALLOYS 

Introduction 

The ability to solidify two -phase aligned structur* s : r i -idc r ?r. ; 

of compositions near a eutectic composition Las importance in t he full 

utilization of this class of composite materials. The early vcrk of 

& Flemings^ claimed that the range of compositions for vhich plane front 

composite growth was possible was restricted not only by constirut ional 

supercooling requirements but also by the requirement t’o:C convection be 

avoided. In fact, their experimental work using vertical upward solid if i^a: ion 

was limited to the Sn-rich side of the Pb-Sn eutectic where the fluid 

density, influenced by temperature and composition, always decreases with 

distance from the interface (height). On the Pb-rich side of the eutectic, 

the concentration gradient in the liquid can cause the density to increase 

with distance near the interface and possibly cause convection (here termed 

[*> 3] 

solutal convection). Subsequently, it was argued that as long as a 

stagnant boundary layer exists whose thickness is large compared to the 

eutectic spacing, plane front growth of the composite should be possible in 

[2 4 5] 

the preseuce of convection. This was in fact verified by experiments ’ ’ 
which demonstrated the growth of Pb-rich Pb-Sn composites. However, Verhoeven, 
Kingery and Hofer^ subsequently showed that significant macrosegregation 
^1C w/v! down the length of solidified samples occurs in samples only 1 w/o 
on the Pb-rich side of the Pb-Sn eutectic, and attributed this macrosegregation 
to solutal convection. 

The purpose of this research is to describe the microstructure and 
macro segregation produced during upward solidification of off-eutectic 
alloys 2-5 w/o from the eutectic composition when the liquid density near 
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the interface increases with distance from the interface and to determine if 


a magnetic field or solidification downward can eliminate solutal convection 
and hence macro segregation. 

Experimental Procedures 

Alloys were prepared from 99.99Z pure Sn and Pb in evacuated and sealed 
borosilicate crucibles. The directional solidification furnace is constructed 
of non-magnt t ic materials and has a bifilar wound cylindrical resistance 
heating element with upper and lower sections and a circulating water cooled 
chill zone. A 3 mm l.D. silica tube containing the desired alloy is pulled 
down through the chill, sealed by O-rings, at constant velocity to accomplish 
directional solidification. The cooling water comes directly in contact 
with the silica tube. Using two separate A. C. proportional temperature 
controllers for the furnace, and water cooled by a heat exchanger, temperatures 
are constant '_o 0.1 °C over an hour and 1 °C over several days. The vertical 
drive is powered by a synchronous motor with a multi-speed gear box. Speeds 
from 0.5 un/s to 0.5 mm/s are possible. The tube is attached to the pull 
shaft using a device which permits the sample to be suddenly pulled down 
about 3 cm to quench the liquid-solid interface. 

The value of the liquid temperature gradient measured for several 
typical experiments with 0.25 mm diameter sheathed thermocouples was 360 + 

10 K/cm. This value was assumed for all other experiments. The furnace is 
designed so that a permanent magnet with a pole diameter of 8 cm and a pole 
separation of 11 cm could be slid around the furnace to apply a 0.1T (1 
kilogauss) horizontal magnetic field. Alternately, an 11 cm diameter solenoid 
could be placed over the furnace-chill to apply up to 0.1T vertical magnetic 
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field. Additionally the furnace could be inverted fot solidification downward. 
In this case a fibrous ceramic plug was used to maintain u small upward 
pressure on the inverted melt to prevent separation of the liquid from the 
freezing solid. 

A specific directional solidification experiment was conducted in the 
following manner. A tube filled with the desired alloy was inserted through 
the chill, with the alloy completely in the chill section. The tube was 
evacuated and the furnace temperature was equilibrated. The sample was then 
driven quickly into the furnace melting about 12 cm of alloy. The tube was 
backfilled with argon and temperatures were allowed to stabilize for about 1 
minute. Excessive holding times at this stage are undesirable because a 
layer of Pb-rich solid solution will slowly solidify and disrupt the uniformity 
of the initial liquid composition. The tube was then pulled at the desired 
rate until about 6 cm had solidified. The remainder of the sample was then 
quenched, delineating the liquid-solid interface in the microstructure. The 
silica was removed with HF acid and the specimens were mounted in a cool- 
setting epoxy for characterization. 

Samples were routinely examined in a longitudinal section (parallel to 
the solidification direction) using standard metallographic techniques. For 
each experiment the shape of the liquid-solid interface was examined to be 
certain that heat flow was unidirectional. In the present experiments, the 
macroscopic liquid-solid interface of two-phase composite growth, when 
achieved, was convex towards the liquid by only 10-20 urn across a 3 mm 
diameter sample. 

Measurement of macro segregation (average composition of the two-phase 
solid) along the length of solidified samples was performed using two methods. 
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x-ray fluorescence and a titration method. The former method, which was 
described previously^ , uses incident monochromatic (22.1 keV) x-rays 
collimated to 1/2 mm x 1/2 ad. The use of x-rays permits deeper penetration 
into the samples and hence better averaging over the two-phase structure 
than does the use of electron probe methods. This data ms used to non- 
destructive^ evaluate the extent of macrosegregation and to indicate those 

[ 8 ] 

samples for which the titration method would be used. The titration method 1 J 
was performed on approximately 2 mm discs cut transverse to the growth 
direction. Accuracy of the titration method is approximately + 0.25 w/o Pb. 

Experimental Results 

Results on the microstructure and macrosegregation of Pb-rich off- 
eutectic alloys solidified upward with and without magnetic field and 
downward are presented here and will be discussed in the next section. 
Directional solidification experiments were conducted vertically upward 
(liquid above the solid) on the Pb-rich side of the eutectic. In all cases, 
solidification begins with the dendritic growth of the Pb-rich phase even 
when the ratio of the temperature gradient in the liquid to the growth 
velocity V substantially exceeds the requirement of constitutional super- 
cooling^ for the bulk composition of the alloy. The volume fraction of 
dendrites decreases slowly as a function of distance as shown, for example, 
in Fig. 1, until the structure is a two-phase composite. The region of 
dendritic growth in most instances is very long compared to D/V (where D is 
the liquid interdiffusion coefficient) and does not correspond to an initial 
transient of the type described by Mollard and Flemings^. This long 
region of dendritic growth will be shown to be a manifestation of the 
macrosegregation which occurs in these alloys due to solutal convection. 
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Table I is a summary of these experiments showing the original composition 
(C q ), growth rate (V), G^/V, D/V, the initial length of the liquid alloy 
(L), the length of the dendritic zone (Z^), and the average solid compositions 
at three positions in the solidified samples: £^(0.0), C^fO.l) and cj* at 

0.0, 0.1 fraction solidified and at the transition from dendritic to composite 
structure, respectively. Fraction solidified, f s , is the ratio of distance 
solidified to the initial length of molten alloy. Other experiments conducted 
with an imposed magnetic field and with solidification downward are included 
in Table 1. 

Fig. 2 shows the measured composition as a function of fraction solidified 
for five 43.0 w/o Pb alloys solidified at 0.5 um/s. [Ihe eutectic composition, c , 
is 38.1 w/o Pb.] Two samples were solidified upward without magnetic field 
with initial melt lengths of 5.2 cm and 11.8 cm resnectively. The initial 
solid composition is very high in Pb content 72 w/o Pb (the boundary of 
the two phase field occurs at 81 w/o Pb) , and falls to the vicinity of the 
eutectic composition as solidification proceeds. The fractions solidified 
f at which the dendritic structure disappears in the two samples are 
essentially identical. Beyond this position the samples continue to solidify 
nondendritically as a two phase composite with changing composition until at 
a f ^ 0.22 the rate of change of composition with distance decreases 
sharply. Because of the excellent agreement of the macrosegregation of 
these two samples as well as the dendrite disappearance position, it is 
apparent that fraction solidified is the proper variable to describe changes 
in composition. This strongly suggests that the dendritic zone is not 
caused by an initial transient but by macrosegregation. Henceforth, we will 
use fraction solidified as the length parameter to describe macrosegregation. 
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Also shown in Fig. 2 are two samples solidified upward with imposed 
horizontal and vertical magnetic fields of 0. IT (1 kilogauss) . The macro- 
segregation and position of dendrite disappearance for these experiments are 
similar to the samples solidified without magnetic field. The final sample 
shown in Fig. 2 was solidified vertically downward (liquid below solid). 

Here a significant reduction in the macrosegregation was observed. Except 
for an initial transient (length ^ 2 D/V) the composition is relatively 
constant. [The slight decrease in composition is caused by thermal diffusion.] 
Furthermore, the structure is completely two phase composite. There is no 
initial region of dendritic growth. 

The effect of a faster interface velocity (2.6 um/s) on macro segregation 
for 43.0 w/o Pb alloys is shown in Fig. 3. Samples solidified upward with 
or without a magnetic field show an initial Pb content of ^ 52 w/o and 
freeze dendritically up to fraction solid of about 0.25. Their macro- 
segregation profiles are slightly different at large f g but certainly 
macrosegregation is quite evident in both. On the contrary the sample shown 
in Fig. 3 which was solidified vertically downward shows very little macro- 
segregation and is two phase composite over most of its length. As shown in 
Table I, the dendritic region in this sample has not been completely 
eliminated by inverted solidification as did occur at 0.5 um/s. This is 
probably due to the closeness of the value of G. /V compared to the critical 

la 

value for plane front growth. (See discussion below and Fig. 7). 

With an initial composition closer to the eutectic composition (40.7 
w/o Pb) macro segregation during upward solidification is reduced compared 
to 43.0 w/o Pb alloys but is still significant as shown in Fig. 4. The 
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sample solidified dendritically to a fraction solid of 0.06. On the 
contrary, the sample solidified downward has a reduced level of macro- 
segregation and is two-phase composite over its entire length. 

Of all the Pb-rich samples solidified upward only one was composite 
over its entire length (with the exception of a 1 nan initial transient). 

The composition of this sample was close to the eutectic composition and it 
was solidified relatively rapidly (40.7 w/o Pb, 5.1 i-m/s) . The macro- 
segregation data for this sample is shown in Fig. 5. 

Discussion of Results 

Upward Solidification - No Field 

The extensive macrosegregation shown in Figs. 1 through 5 occurs only 

for alloys on the Pb-rich side of the eutectic when freezing vertically 

upward. Experiments on the Sn-rich side of the eutectic under conditions 

identical to those used in the Pb-rich experiments exhibit only a relatively 

small amount of macro segregation which can be quantitatively accounted for 

by thermal diffusion alone. An example of measured macrosegregati. i on the 

Sn-rich side is shown in Fig. 6. The data agrees with previous work of 

rn] 

Verhoeven, Warner & Gibson , and can be fit using their equation (17) with 

2 

a thermal diffusion coefficient, D, of 2.1x10 cm /(s*K), a diffusion co- 

-5 2 

efficient of 0.6x10 cm /s and no mixing near the liquid-solid interface (6 
large). This absence of significant macrosegregation on the Sn-rich side in 
experiments conducted under conditions identical to those on the Pb-rich 
side is strong evidence that solutal convection (and not thermal convection 
caused by radial temperature gradients) is the cause of the extensive 
macrosegregation on the Pb-rich side of the eutectic. 
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It is useful to describe a qualitative model for the macro segregation 
and resultant changes in microstructure for Pb-rich alloys. Solutal 
convection is sufficiently vigorous to disrupt the normal diffusion con- 
trolled (Sn-rich) solute field near the liquid-solid interface regardless of 
whether it <s planar or dendritic. This solute is transported by the 

convection to the bulk liquid far from the interface which is almost always 

[9] 

well mixed* by the radial temperature gradients which exist there . 

Hence the bulk liquid, which initially was of composition C q on the Pb-rich 
side of the eutectic, gradually becomes more Sn-rich. This causes a slow 
change in the solid to compositions which are also more Sn-rich as solidi- 
fication proceeds. The transition from dendritic to composite structure 
appears to occur when the solid composition is close enough to the eutectic 
composition that plane fro:.t composite growth becomes stable as documented 
below. After plane front solidification is established, macrosegregation 
continues until the bulk liquid composition approaches the eutectic 
composition. Because the average composition in the liquid at the interface 
must be very close to the eutectic composition for two-phase plane front 
growth, the average composition gradient in the liquid vanishes when the 
bulk liquid composition reaches the eutectic composition. At this point 
solutal convection presumably ceases. This situation is reflected by the 
sharp break at f 'v 0.22 of the macrosegregation curve in Fig. 2. The fact 
that the solid composition at this point is ^ 44 w/o Pb and not equal to the 


* In the present experiments the composition of the liquid, determined 
after quenching, was constant to within + 0.25 w/o Pb. The large 
sample size 2 mm} used for chemical analysis prevents detection of 
the solute boundary layer. 
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liquid composition (3ft . 1 w/o Pb) is caused by a thermal diffusion effect^. 
From tq. (15) of Ref. 9, the average composition gradient vanishes when 
the average solid composition, C g , is C £ (1+D'G L /V) , which for this alloy 
and these growth conditions is A3. 9 w/o Pb. Macrosegregation after 
this point is greatly reduced and is caused only by thermal diffurion. 

A quantitative description of the location of dendrite disappearance is 
possible. Fig. 7 is a plot of G T /V vs. the composition of the solid at the 
position where the structure changed from dendritic to composite. It 
resembles the expected form of the relation between G /V vs. composition for 
the stability of plane front solidification, i.e., it falls slightly below 
the requirement of the simple constitutional supercooling criterion 

G t /V = - m(C -C )/D 
Li t s 

where m = 3.2 K/(w/o) and is the liquidus slope of the Pb-ricl. phase. It 

[ 2 ] 

also agrees with the data of Davis & Fryzuk on zone melted off-eutectics 
where macrosegregation is avoided by the fact that solid of constant com- 
position is fed into the melted zone. Hence the transition from dendritic 
to composite structure occurs when macrosegregation has caused the composition 
of the solid to change to a level which satisfies a stability criterion for 
plane front growth of a composite. 

Quantitative description of the observed macrosegregat ion is more 

f 6 91 

difficult. Verhoeven et al. 1 ’ , using a thermotransport-corrected Burton- 

Prim-Slichter type (stagnant film) analysis of plane front composite 
growth, have formulated an expression to predict macro segregation. The AO. 7 
w/o Pb sample run at 5.1 um/s is dendritic for an extremely short distance 
(long compared to D/V) and hence this analysis is applicable. The pre- 
viously given value of the thermal diffusion coefficient D % and a value of 
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90 pm for 6 (the thickness of the unmixed zone near the interface) gave the 

best fit to the measured data In close agreement with 6 values of Verhoeven 
[9] 

et al. . This calculated curve of the nacrosegregation is shown super- 
imposed on the data of Fig. 5. In this experiment, it was determined 
experimentally that the liquid composition remained on the Pb-rich side of 
the experiment up to the quench po'at and solutal convection probably existed 
throughout this experiment. The calculated value of the bulk liquid composition 
using this value for 5 also agreed to within 0.25 w/o with the measured 
composition of the quenched liquid. The success of this approach during 
composite growth is somewhat suprising given the fact that the driving force 
for solutal convection only exists in the diffusion boundary layer. On the 
contrary, samples With iuug regions -f dendritic growth cannot be fit with 
this equation In the data shown in Fig. 3 a substantial number of data 
points exist in the dendritic zone of the sample and a marked change in 
slope occurs as tne volume fiaefio.. tz dendrites decreases. Clearly the 
existence of the dendritic interface influences both the driving force for 
solutal convection and the geometry for the flow. Fluid flow during 
dendritic solidif ication^^ is known to require a different approach to 
macrosegregation than a stagnant film approach. Quantitative description of 
the macrosegregation curves during dendritic growth is a topic for future 
research. 

It is significant to note several aspects of solutal convection observed 
in this research. Fig. 8 shows an optical micrograph of a quenched liquid 
solid interface of a 43.0 w/o Pb alloy freezing at 2.6 pm/s. The sample was 
quenched at f = 0.1, so solutal convection was present at the time of 
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quenching. The Pb dendrites are uniformly spaced across the sample an i the 
eutectic interface appears quite flat. It seems unclear that a single (or 
double) convection cell would lead to such a microstructure. 

Additional insight into the nature of solutal convection is seen in 
Fig. 2, where the sharp break in the data at f c 'v 0.22 marks the transition 
from macrosegregation caused by solutal convection to macro segregation 
caused by thermal diffusion. This sharp transition from regimes of con- 
vection to no convection may be an indication of the extremely low concentration 
differences (i.e., deviation from eutectic composition) necessary for the 

onset of convection. Low concentrations for the onset of solutal convection 

[ 12 ] 

during single phase solidification have been predicted theoretically 

For example, for vertical solidification of single phase Pb containing Sn, 

at an interface velocity of 1 pm/s and a temperature gradient of 200 K/cm, 

—4 

solutal convection will occur at all compositions above 3x10 w/o Sn. 

Another suprising feature of solutal convection is the failure of the 

small (3 mm) diameter container to suppress convection. The results of 
f 121 

Coriell et al. indicate that a solute Rayleigh number based on a length 

of D/V is significant for the prediction of the onset of solutal convection 

during solidification. The fact that D/V is smaller than the radius of the 

tubes used in this study may explain this observation. 

The presence of growth rate fluctuations has been observed in some 

cases in the present study. However, as can he seen for example in Fig. 1, 

many samples have been obtained vhich contain no fluctuations. Hence, it is 

our belief that these fluctuations are caused by artifacts, such as for 

[13] 

example, by nonsteady thermal contractions . Solutal convection seems to 
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be sufficiently steady so as to not cause growth rate fluctuations. 

Upward Solidification - Magnetic Field 

Static magnetic fields have long been known to influence the flow of 
electrically conducting liquids . Chandrasekhar calculated the 
influence of a vert:’ al magnetic field on suppression of thermal convection 
in an infinite horizontal layer of fluid heated from below. Other theo- 
retical work has been performed for different geometries. Experiments on 
liquid mercury confirmed the inhibition of thermal convection by a 
vertical magnetic field. During solidif icatio n , liquid temperature 
fluctuations and associated growth rate fluctuations and solute banding 
caused by time- dependent convection during the horizontal growth Cj. crystals 
were eliminated by imposition cf a vertical magnetic f ield^’^’^ . 
However, time independent laminar flow may still have been present. Because 
of this previous work, attempts were made in the present study to suppress 
solutal convection and the concomitant macrosegration by the imposition of 
a static magnetic field. 

For most experiments a vertical field of 0.3T (1 kilogauss) was used 
based on the following analysis from Chandrasekhar . The fluid elements 
are subject to the Lorentz force which is proportional to curl H x H, where 
H is the magnetic field in the fluid. Clearly, if H is constant the Lorentz 
force vanishes, and the magnetic field does not affect the fluid motion. 

Even though the applied magnetic field is constant, the magnetic field in 
the fluid is not in general constant. In fact, a constant H will satisfy 
the basic partial differential equation for H if and only if curl (u x H) 
van Lshes, where u is the fluid velocity. Clearly if u is parallel to H, 
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then u x H vanishes, and the magnetic field has no effect on the fluid 
motion, i.e., if the fluid motion does not cross the magnetic field, there 
is no effect. More generally, however, a constant H is possible if and only 
if (H*V)u = 0 for an incompressible fluid, which reduces to (?u/3x.) = 0 
for a magnetic field only in the di*ection. Thus, a constant magnetic field 
and a vanishing Lorentz force will occur if the fluid velocity is independent 
of the coordinate in the direction of the magnetic field. Such a case 
occurs for flow between two infinite horizontal plates heated from below 
where the critical Rayleigh number is unaffected by a horizontal magnetic 
f xeld. 

Experimental results shown in Figs. 2-3, show little or no effect of a 

0.1T (1 kilogauss) magnetic field on reducing the macrosegregation produced 

by solutal convection. As was described above, the solutal convection 

present in these experiments is sufficiently steady that growth rate 

fluctuations do not occur. Hence the nature of the convection present in 

these experiments ijs quite different from that observed by Utech & Flemings 

during horizontal solidification. Recent theoretical work may explain the 

[ 12 ] 

ineffectiveness of a magnetic field on damping solutal convection . For 
example, for vertical solidification of single phase Pb containing Sn, at an 
interface velocity of 1 ;jm/s, a temperature gradient c 200 K/cm and a 
imposed vertical magnetic field of IT (10 kilogauss), solutal convection 

_3 

will occur at all compositions above 3 x 10 w/o Sn. 

Downward Solidification 

As shown by the data in Figs. 2-4, solidification of Pb-rich off- 
eutectic alleys downward greatly reuuces macrosegregation. In fact the 
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remaining macro segregation in these experiments can be quantitatively 
accounted for by thermal diffusion. In this geometry the Sn-rich solute 
layer near the interface is above the Pb-rich liquid away from the interface 
and for the conditions used in these experiments, creates a stable density 
gradient near the interface even though the temperature gradient is destabilizing. 
Outside the solute diffusion layer the density gradient is destabilizing. 

At first glance the reason for the reduction of macrosegregation during 
downward solidification appears obvious. Solutal convection is eliminated 
because of the stable density gradient near the interface and hence diffusion 
controlled solidification is not interrupted. However, the large vertical 
temperature gradient present during these experiments would tend to promote 
extensive thermal convection farther from the interface. The results indicate, 
however, that such flow if it exists, is sufficiently distant from the 
diffusion boundary layer as to be ineffective in causing significant macro- 
segregation. 

A second rationale could be that the small tube diameter (3 mm) is 

effective in greatly reducting the thermal convection far from the interface. 

One can calculate a thermal Rayleigh number R - a gG r^/v< for this geometry, 

where a is the expansion coefficient, g the acceleration of gravity, the 

temperature gradient in the liquid, r the tube radius, v the kinematic 

[121 

viscosity and < the thermal diffus’vity. For typical parameters of the present 
experiments R = 65, which is very close to typical critical Rayleigh numbers 
for closed containers as summarized by Verhoeven^*^ . for larger tubes 
thermal convection away from the interface may become more effective in 
causing macro segregation and hence it is anticipated that suppression of 
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macro segregation by solidification downward may be very sensitive to tube 
diameters. This effect was not, however, investigated in the present study. 

Conclusions 

1. Solute gradient induced convection can cause extensive macro- 
segregation in off-eutectic alloys. 

2. Macro segregation is frequently sufficiently large to prevent plane 
front solidification of the alloy. 

3. The transition from dendritic to composite structure occurs when 
the composition of the solid is close enough to the eutectic 
composition to satisfy a stability criterion based on G^/V. 

4. Horizontal or vertical magnetic fields of 0.1T (1 kilogauss) do 
not prevent macro segr egation of oilf-eutectic alloys caused by 
solutal convection. 

5. For si^all diameter samples, inverted solidification eliminates 
macrosegregation caused by solutal convection. 
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B. CALCULATION OF CONDITIONS FOR THE ONSET OF SOLUTAL CONVECTION 
Introduction 


During the vertical directional solidification of a binary alloy at 
constant velocity V, the steady state concentration field is exponential 
in form, i.e., c =* c^+Cc^/k) (1-k) exp (-Vz/D), where c is the concentration 
of solute in the liquid at a distance z from the solid- liquid interface, 
the distribution coefficient k is the ratio of the solute concentration 
in the solid at the solid-liquid interface to the solute concentration in 
the liquid at the solid-liquid interface, c b is the solute concentration 
as z -*- 00 , and D is the liquid diffusion coefficient. The temperature 
field is also exponential but with decay length k/V, where < is the thermal 
diffusivity of the liquid. In general, k » D so that the temperature 
gradient G L in the liquid is essentially constant in the region in which 
the concentration gradient is significant. It is well known that as the 
concentration is increased for fixed V and the planar solid-liquid 
interface becomes morphologically unstable. The constitutional supercooling 
principle predicts instability when mG > G , where m is the slope of the 

C L 

liquidus line and G = (c V/Dk) (k-1) is the concentration gradient at the 

c °° 

planar solid-liquid interface. In addition to morphological instability, 
there is the possibility of a convective instability. For solidification 
vertically upwards (positive z-direction) the thermal gradient alone 
results in a negative density gradient (stabilizing) while the solute 
gradient alone can cause a negative (stabilizing) or positive (destabilizing) 
density gradient depending on the density variation with solute concentr~‘ * on 
and the distribution coefficient. When the solute gradient gives rises 
to a positive density gradient, e.g., a solute which is rejected at the 
solid-liquid interface and which is lighter than the solvent, a convective 


i 

5 
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instability will occur if the solute concentration is sufficiently large 

that the destabilizing solute field overcomes the stabilizing influences of 

the thermal field and viscous effects. Thus, for fixed values of velocity V 

and temperature gradient G , as the solute concentration increases, one 

expects an instability to occur, i.e., for suff iciently small solute 

concentrations, the planar interface is stable and the horizontal components 

of the fluid flow velocity vanish (a vertical flow due to the volume change 

on solidification is always present), but for solute concentrations above 

some critical value the interface and the fluid flow field are unstable. 

[12 21 22 ] 

In our previous publications ' ’ , a linear stability analysis and 

numerical algorithms have been described which allow the calculation of the 
critical concentration for the onset of convective and morphological in- 
stabilities during the directional solidification of binary alloys. In 
addition, a number of results were obtained for the specific case of the 
solidification of lead containing tin. In this report we combine the 
previous results and new results with particular emphasis on small constant 
gravitational accelerations. 

Results and Discussion 

He recall that the dependence of the perturbed temperature, con- 
centration, fluid velocity, and interface shape on the horizontal coordinates 

x and y and the time t is of the form exp [ot + i(u> x+u y)], where in and co 

x y x y 

are spatial frequencies. If the real part, o^, of the time constant a is 

positive for any values of and the system is unstable. If < 0 

for all values of u and w , the system is stable. Ue refer to instabilities 

x y’ 
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corresponding to large values of u (small wavelengths, A « 2ir/ui) as 
morphological instabilities and to instabilities corresponding to small 
values of u (large wavelengths) as convective instabilities. 


The basic results for the critical concentration c** above which 

■ o> 

instability occurs as a function of growth velocity are shown in Fig. 9. 

We have previously discussed the curves corresponding to gravitational 

2 -4 -6 

accelerations of g £ = 980 cm/s and 10 g g . For g * 10 g g , the con- 
vective instability occurs at a lower concentration than the morphological 
instability only if V < 1.3 pm/s. The concentration above which the 

2 2 k 

interface is unstable as a function of the spatial frequency w * ( / 

of the perturbation is shown in Figs. 10 and 11 for V * 1.2 and 1.3 pm/s, 

respectively. The critical concentration c£* corresponds to the m-tn-timnn 

value of the concentration as a function of u. The solid curves on the 

right of Figs. 10 and 11 correspond to morphological Instability while the 

solid curves on the left correspond to convective instability. Along the 

curves the real part of a vanishes; along the solid curves o vanishes 

while along the dashed curve, which connects the convective and morphological 

branches, the imaginary part o^ of a is non-zero. The value of is 

shown in the insets in Figs. 10 and 11 and attains a maximum value of less 
-4 -1 

than 2.0(10 )s . For V * 1.2 pm/s, the minimum value of c b is 4.61 wt. 

Z at u * 1.35 cm ^ or a wavelength A = 4. 7 cm. There is a relative minimum 
on the morphological branch at u * 204 cm ^ with c^ = 6.47 wt. Z (not 
shown in the figure). For V * 1.3 pm/s, the minimum value of c^ is 5.98 
wt. Z at u = 210 cm ^ (not shown in the figure) and a relative minimum 


with c 


7.15 wt. Z at u * 1.4 cm 


-1 
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In our previous publications, it was shown that the critical concen- 
trations for convective instability could be correlated in terms of the 
Rayleigh number R' and the solute Rayleigh number S' based on the length 


(D/V), i.e., 

R' - -gaG^D/V) 4 /^ 

S' - -go G (D/V) 4 /vD - (ga /vDk)(l-k)(D/V) 3 c , 
c c c 00 

where a is the thermal coefficient of expansion, v is the kinematic 

viscosity and a is the solutal coefficient of expansion. To fit our 

previous results we used an expression of the form 

a 5 

S' + a^R') = a 2 


( 1 ) 


or in dimensional form 

c a = a 3 V 3 /g + a 4 (V 3 /g)(gG L /V A ) 5 , (2) 

where the a. are constants independent of G T , c , V and g. The a. were 
determined by minimizing the differences between the numerical values of 
c a and the values obtained from the above expression. Using the previously 
determined values of the a we calculate critical concentrations of 
0.32, 2.48, 4.27, and 5.42 wt. Z for V * 0.5, 1.0, 1.2, and 1.3 pm/s, 
respectively, for g * 10 ^ g g ; these are in reasonable agreement with the 
numerically calculated values of 0.31, 2.39, 4.61, and 7.15 wt. Z for V = 
0.5, 1.0, 1.2, and 1.3 pm/s. 

In our previous calculations for g * g^, we obtained oscillatory in- 
stabilities with 'p 0.1 s 3 or a period of the order of 60 s. For 

-6 -4 -1 4 

g * 10 g p ^ 10 s or a period of the order of 6(10 )s. This 

suggests that the fluid motions are very slow at reduced gravity levels. 
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To further verify this, we carried out a few calculations with ^ 0. 

For V ■ 1.2 pm/s and g * 10 ^ g g , a concentration of 7.68 wt. Z is required 
-4 -1 

for o^-lO s on the convective branch. Thus, although o > 0 for 

c d > 4.61 wt. Z for V * 1.2 un/s, the anplitude of the convective in- 

4 

stability for c^ » 7.68 wt. Z increases very slowly, requiring 10 s to 
increase by a factor of e. In contrast, the morphological instabilities 
are only slightly affected by reduced gravity. For example, for V - 1.2 

—1 —4 

pm/s and u) » 205 cm , the concentrations corresponding to * 0, 10 , 

—2 —1 —*6 
and 10 s are 6.47, 6.50 and 11.0 wt. Z, respectively, for g - 10 g £ 

-4 -2 -1 

whereas for g = g g , * 0, 10 and 10 s correspond to concentrations 
of 6.47, 6.49 and 11.2 wt. Z, respectively. We emphasize that these 
calculations for g = .10 ^ assume that g is constant, and since the 
convective instabilities develop very slowly, it seems unlikely that this 
is a valid assumption for space experiments. 

Liquid metals have low Prandtl numbers, e.g., the Prandtl number 
Pr = v/< = 0.0225 for lead. In contrast, a large class of liquids have 
much larger Prandtl numbers, e.g., Pr ^ 10 for water. In order to study 
the effect of Prandtl number, we have carried out some calculations using 
the properties of lead-tin alloys except that the thermal conductivities 
and thermal dif fusivities of both liquid and solid were 100 times lower 
than the actual conductivities and dif fusivities. Thus, the Prandtl number 
for this hypothetical alloy is 2.25; for fixed and V the thermal Rayleigh 
number is also increased by a factor of 100. The tin concentrations at the 
onset of convective instability as a function of growth velocity for * 

200 K/cm for lead-tin alloys (Pr » 0.0225) and the hypothetical alloy 
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(Pr * 2.25) are shown in Fig. 12. At low velocities the critical con- 
centrations are a factor of 50 greater for the liquid with reduced thermal 
properties than for the lead alloy. At higher velocities, the two curves 
approach each other and the thermal properties have little effect on the 
critical concentration. This behavior can be understood qualitatively 
from eq. (1) through the effect of the reduced thermal properties on the 
Rayleigh number; in general, however, we expect the a^ to depend on the 
Prandtl number. 

A few additional calculations of the effect of magnetic fields and 
thermal gradients were carried out. These together with previous results 
are presented in Figs. 13 and 14. For discussion of these results, see the 
previous publications. 

The present linear stability analysis and numerical algorithm allows 
the determination of the onset of convective and morphological instability 
during directional solidification. We plan additional calculations for 
conditions and materials relevant to experimental investigations. For a 
particular double diffusive system with linear temperature and solute 
gradients, Baines and Gill^^ have shown that the spatial frequency u, 
corresponding to the most unstable mode, increases very rapidly as the 
concentration exceeds the critical concentration for the onset of in- 
stability. Calculations of o as a function of w for fixed concentration 
c m will allow us to determine whether similar behavior occurs during 
directional solidification. Mathematical analysis and algorithm development 
for fluid flow occurring in the non-linear regime (after the onset of 
instability), with particular emphasis on solute distribution, will be 
initiated in the near future. 
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Table I. Summary of Experiments 
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Longitudinal section of dire.ctionally solidified sample 
showing the change in vol. me fraction of Pb dendrites (dark) 
as a function of distance (fraction solidified). Solidification 
proceeded from bottom to top. At the extreme top there are 
nc dendrites. Solidification velocity is 1.3 um/s - original 
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Fig. 2. Average solid composition versus fraction solidified for 

43.0 w/o Pb alloys solidified at 0.5 pm/s showing extensive 
macro segr egation during upward sc .idification, negligible 
effect of magnetic field and large effect of inverted solidi- 
fication on reduction of macrosegregation. 
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Fig. 4. Average solid composition versus fraction solidified for 
40.7 w/o ?b alloys solidified at 1.3 ym/s. 
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upward at 5.1 pm/s. The curve is the best fit of the data 
using only 6 as a parameter in equation (17) of Ref. (9). 
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Fig. 7. A plot of G. /V versus solid composition at the position of 
dendrite disappearance for all of the samples. The data 
follows the general trend of plane front interface stability 
for off-eutectic alloys. 





dP/di -C 


V /x m /$ 


The critical concentration c** of tin above which instability 

O) 

occurs as a function of the velocity v of directional solidi- 
fication of lead-tin alloys for a temperature gradient C in 
the liquid of 200 K/cm. The solid curves represent gravita- 
tional accelerations of g = 980 cm/s ,10 g , and 10 g . 
The upper line with negative slope represents the onset of 6 
morphological instabilities; the nearly parallel dashed-dot 
line labeled Oi>/3z) * 0 represents the neutral density 
criterion. The dotted extension of the curve labeled g 
corresponds to oscillatory instabilities. e 
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Fig. 10. The concentration at the onset of Instability during direc- 
tional solidification at V = 1.2 ym/s of a dilute alloy of tin 
in lead as a function of the spatial frequency u of a sinusoidal 
perturbation; the temperature gradient G T in the liquid is 200 
K/ci^and the constant gravitational acceleration is 980 (10 - ) 
cm/s . The solid curves mark the onset of non-oscillatory 
instabilities (o^=0); whereas the dashed curves mark the onset 
of oscillatory instabilities (the value of a is given in the 
inset) . 
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Fig. 11. The concentration at the onset of instability during direc- 
tional solidification at V = 1.3 yra/s of a dilute alloy of tin 
in lead as a function of the spatial frequency u of a sinusoidal 
perturbation; the temperature gradient C. in the liquid is 
K/cm^and the constant gravitationax acceleration is 980 (10 ) 

cm/s . The solid curves mark the onset of non-oscillatory 
instabilities (a =0); whereas the dashed curves mark the onset 
of oscillatory instabilities (the value of is given in the 
inset) . 
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TIN CONCENTRATION w/o 



I 10 100 

INTERFACE VELOCITY ,i/s 


Fig. 12. The critical concentration of tin above which instability 

occurs as a function of the velocity of directional solidi- 
fication for a temperature gradient in the liquid of 200 K/cm. 

The curve with Prandtl number Pr = 0. 0225 corresponds to a 
lead-tin alloy whereas the other curve corresponds to calculations 
for a hypothetical lead-tin alloy with thermal properties 
reduced by a factor of 100, i.e. , Pr = 2.25. 
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i Fig. 13. The critical concentration of tin in lead above which in- 

; stability occurr as a function of the temperature gradient in 

< the liquid for solidification velocities of 1.0 and 20.0 pm/s. 
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Task 3 


A Thermochemical Study of Corrosive Reactions in Oxide Materials 

H. S. Parker 

Ceramics, Glass and Solid State Science Division 
Center for Materials Science 

Summary 

Phase assemblages of relerted compositions on the joins 
K[Fe jSi j]02”Si02 and KFe02 HO 2 have been determined after heat 
treatment at oxygen partial pressures of 10~^ and 10 -1 ^ atm and 
temperatures in the 800°-1400°C range. Excessive reaction between 
platinum containers and melts containing some ferrous iron was noted. 
A comparison of conventional and containerless techniques for phase 
equilibrium studies is presented. 
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Introduction 


The objectlv. a of this task are to Investigate the nature, extent 
and llnliations on experimental studies Imposed by sample-container 
reaction at high temperature and to evaluate containerless techniques 
for Belt and solid state Investigations. Alkali and transition elements 
ere among the most reactive and are of great technological Importance 
because of their application in the fields of MHD, fuel cells, 
catalysts and s^xid state electronic and ionic conductors. Reaction 
with the con r can severely compromise the bulk composition of the 
sample and adversely affect the electrical and chemical properties. 

TMs, together with the Inability in general to electromagnetically 
levitate ceramic materials, makes low gravity experiments attractive 
for containerless study of mch materials. 

The general approach la to select a system of current technological 
importance and Investigate: a) the maximum temperature limits for the 
container material and characterize the reaction products, b) determine 
the equilibrium phase assemblages in selected portions of the systems 
within the limitations Imposed by a) , aud c) in collaboration with MSFC 
provide specimens for levitation experiments at MSFC and characterize 
“•nd correlate results of levitation experiments with container experiments 
to evaluate problems associated with lov gravity experiments. 

The system being Investigate! is the O-iron oxide-SiO^ system. 

This syxtf-' *s of aujor Interest to the magnetohydrodynamics (MHD) 
effort, where these oxides are present as components of seed and coal slag. 
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The fleet portion of this system Investigated daring the course of this 
project was compositions on the KFe02*‘ n Fe20^ n binary. Due to the high 
temperatures involved (T > 1500* C) platlma was chosen as the container 
material. As a result of this study. It was reported [1] that, contrary 
to the expected result, platinum was a suitable container at temperatures 
as high as 1600*C, provided most of the Pe was kept in the 34- state. 

Also, the existence of a 8" '-alumina structure type phase was reported. 

As an extension of this study, selected compositions In the K^O- 
lron oxlde-SlOj system have been examined at lower partial pressures of 
oxygen, where divalent Iron Is present. 

Experimental Procedures 

The preparation of compositions for study and characterization methods 
following heat treatment has been described [1,2]. The same methods were 
utilized for this portion of the study, except that all calcinations and 
beat treatments were performed in controlled partial pressures of oxygen. 

In order to provide atmosphere control during heat treatment of 
-eclmens, a gas blending find purifying train was assembled. Mass flow- 
meters were used to establish desired nixing ratios of Inert gas (argon), 
95 N 2 / 5 H 2 , N 2 and 0^. The Inert gas was further purified by passing 
over hot titanium sponge at approximately 800*C when desired. The actual 
oxygen partial pressure In the furnace atmosphere was continuously 
monitored using a stabilized zlrconla oxygen sensor. The furnace chambers 
consisted of impermeable high purity alumina tubes provided with suitable 
metal end fittings for gas flow and temperature measurement, with silicone 
rubber gaskets. 
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Because of the lower temperatures involved , the choice of container 
materials included iron, gold and platlnun. The use of iron as a con- 
tainer aatei.il must be limited, however, to oxygen pressures determined 
by the Fe-FeO equilibrium, of the order of 10*^ to 10 -1 ^ atm at the 
temperatures Involved In this study. Gold containers were utilized 
for heatings at higher oxygen pressures and temperatures up to the 
1000 > -1050*C range. For heatings above this temperature, platinum 
containers were used. 

Experimental Results 

The compositions, temperatures and oxygen partial pressures investi- 
gated are shown in Table 1. Selected initial compositions on two joins 
were examined, the K(Fe 5 Si ^) 2 ~Si0 2 join in the K 2 0-"Fe0 ,, -Si0 2 system 
and the analogous KFeOj-SiOj join in the K^O-Fe^^-SiC^ system. Thus, 
all final compositions fall between these bounding planes in the system 
K-Fe-Si-O. The compositions in the K(Fe ^Si 5 )0 2 ~Si0 2 join are an 
extension to lower silica content of the region studied earlier by 
Roedder [3]. The KFe0 2 -S10 2 join has been studied in air by McDaniel [4]. 

When heated at temperatures ox 800*C and o. .- gen partial pressures 
of 10 atm, all compositions on the K[Fe 5 Si 5 ]0 2 -Si0 2 join yielded a 
phase similar in structure to the high temperature orthorhombic form of 
KFeSiO^ previously reported [4,5], In air, this orthorhombic KFeSiO^ 
transform, eversibly tc a low temperature hexagonal form at about 945°C. 
Apparently, the reduced oxygen pressure during these experiments increases 
the stability region of the orthorhombic structure. Reheating in air at 
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800°C of a sample of 2K(Fe ,-Si ^)0 2 P reviou ®ly heated at 800*C In 
—6 

P ■ 10 atm (Table 1) resulted In oxidation and formation of the 
°2 

expected hexagonal , low temperature structure of KFeSiO^. Heating at 
lower partial pressures of oxygen, 10 ^ atm, showed decomposition of 
the orthorhombic structure, with the appearance of free Si0 2 In the 
x-ray diffraction pattern. 

Compositions on the KFeO^-SiOj join, prepared by calcining in air 

to yield trivalent iron, showed generally similar results to those 

initially containing divalent iron in that phases having structures 

similar to hexagonal and orthorhombic KFeSiO^ were observed after heating 

at reduced oxygen pressures. The initial composition AKFeO^tSSiO^ 

is known [4,5] to have a tetragonal structure in air, at temperatures 

above 850°C, however, apparently not stable at partial oxygen pressures 

of 10~* atm. Melting points of compositions on this join are lowered 

100°-150°C by heating at the oxygen pressures shown in Table 1 as 

compared to the same compositions heated in air [4,5]. 

The use of gold containers for low temperature experiments was 

generally satisfactory. In cases where melting occurred, minor gold- 

sample reaction was evident. For temperatures above the gold limit 

where platinum was used, reaction between the molten specimen and 

platinum was pronounced. Attempts to prepare homogenous glasses for 

equilibrium recrystallization experiments at oxygen partial pressures 

_g 

of 10 atm and temperatures in the 1200°-1400°C range were unsuccessful. 

Three compositions were tried: KFe + ^9iO^, 2K[Fe + jSi 5 ]C> 2 :Si0 2 and 

i 

! 



+2 

5K[Fe 5 Si 5 J0 2 :Si0 2 . Reaction between the platinum containers and 
the molten specimens was observed In all three compositions. Consider- 
able embrittlement of the platinum due to iron pickup was evident. 

In the case of the 1400*C heating, traces of the melt were observed 
on the outside of the container suggesting diffusion through the 
container wall. 

Conclusion and Discussion 

Phase assemblages of selected compositions on the joins K[Fe ^Si^JO^-SiC^ 
and KFeO^-SlO^ have been determined after heat treatment of specimens 
at oxygen partial pressures of 10 ® and 10 ~^ atm. The use of gold as a 
container appears to be suitable for lover temperatures where no liquid 
is present. Attempts at higher temperatures to prepare homogenous glasses 
for equilibrium recrystallization experiments were unsuccessful due to 
excessive reaction between the platinum container and feirous iron in the 
melt at these oxygen pressures. 

Discussions with Dr. Oran of MSFC have shown tha v it is not feasible 
to acoustically levitate in a 1-g environment and simultaneously heat 
specimens of these materials to the required temperatures, either for 
preparation of homogenous glasses or for melting point studies. An 
alternate technique is the use of aerodynamic jet levitation and it is 
suggested that this possibility be explored further. 

A comparison of experimental conditions and requirements for con- 
ventional and containerless phase equilibria studies as well as 
anticipated problem areas are given in Tables 2 and 3. 
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Table 2. Comparison of Experimental Conditions for Phase Equilibria 
Studies. 


CONVENTIONAL 


CONTAINERLESS 


Specimen Pre-reacted , powder, packed 

into tube, tube sealed. 

Mass 'V0.1 to 0.3 g 


Maximum Heat- No restrictions, plunge 

ing Rate into hot furnace 


Temperature 

A. Measurement Thermocouple or optical 
pyrometry ±1° to ±5 # C 


3. Control 


±2° to ±5°C 


C. Time at 
Temperature 

1. Subsolidus 

Experiments hrs to days 


Pre-reacted pcvder, formed 
into sphere or right 
cylinder. Mass M3.1 to 
0.3 g 

Limited only by thermal 
stress resistance of 
specimen 


Two color pyrometry or 
non-contacting method 

Better than ±25°C at 
final temperature 


Not feasible at present 


2. Solidus Minutes to 1-2 hrs after 

and isothermal equilibrium 

Liquldus 
Determination 


Minutes to 1 hr after 
isothermal equilibrium 


Cooling Rate ^nch into water 


As rapidly aj possible 


Characterization 
and Evaluation 


In all cases, specimens are examined post-test 
by optical and scanning electron microscopy 
and x- rav ■’if fraction to determine phase 
assemblage, partial or complete melting. 

Where appropriate, chemical analyses of 
specimens are performed. 



Table 3. Some Problem Areas Associated with Containerless Experiments. 


1) Heating technique for 1300° - ld00°C range 

1) radiant wall furance 

2) laser 

3) arc -image furnace 

2) Temperature measurement and control 

1) two-color pyrometer 

2) other non-contacting method 

3) Thermal gradients in specimen must be as small as possible 

• 4 ) Charges in bulk composition of specimen must be measured or 
estimated from vapor pressure measurements 

5) Density of specimen will change with sintering and melting 
during experiment 

6) Electrical conductivity will vary with both temperature and 
composition 

In the case of low gravity experiments, electromagnetic positioning 
is an alternate technique applicable tc a more limited class of oxide 
materials. As part of an investigation in the area of refractory oxides 
having good electronic conductivity, a g-oup of NBS Investigators has 
recently reported [6] conductivities as high as 1 ohm - *cm~* in the 100°- 
250°C range in calcium-doped yttrium chromites. Melting points ate in 
excess of 1750°C, making the choice of container material severely limited 
for high temperature property measurements. The conductivities appear to 
be large enough to make electromagnetic positioning feasible. 
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Task 4 


} 

Thermodynamic Properties of Refractory Materials 
at High Temperatures 

Jack H. Colwell 

Chemical Thermodynamics Division 
Center for Thermodynamics and Molecular Science 

Summary 

The current status of the program to obtain high temperature (1500-4000 K) 

heat capacities of materials while they are levitated in a space environment 

is presented. The heat capacities are to be derived from the simultaneous 

measurement of the radiant heat loss and the rate of temperature change of 

spherical samples as they cool freely. Details of the thermal gradients to 

be expected in both metallic and nonmetallic samples have been examined. In 

metallic samples the gradients, in general, will be tolerably small and cause 

little difficulty. In nonmetallic samples the gradients will be large but 

by making measurements on samples of different size, it should be possible to 

get moderately precise values of botn the heat capacity and the thermal 

conductivity. The upper temperature limit of the experiments wilJ be determined 

by sample vaporization and will occur when the vapor pressure reaches the 

-4 

vicinity of 10 Pa (10 atm) . The program for a ground-based simulation of 
the free-cooling experiment is discussed. 
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Introduction 


The purpose of this research is to measure the heat capacities of 
refractory materials at high temperatures, 1500-4000 K, by using the low- 
gravity environment of space to maintain the samples without physical 
contact. The samples will be spherical, fabricated in that form as solids 
and maintained in that form as liquids by their surface tension. They will 
be heated to a high temperature and then allowed to freely cool by radiation 
in a cold-walled vacuum chamber. 

As the sample cools its temperature will be determined by radiation 
pyrometry and the total radiative heat loss will be determined simultaneously 
using a wide-band pyroelectric detector. The ratio of the heat loss to the 
rate of change of temperature will be a measure of the heat capacity. 

To give an indication of the nature of the experiment, two pertinent 
quantities are presented in Fig. 1. The rate of cooling of a sphere by 
radiation is obtained, assuming no thermal gradients, by equating the rate 
of heat loss 

dH = -47rR 2 £0T 4 dt (1) 

with the thermal capacity 

dH = (4/3)7TR 3 cpdT (2) 

which yields 

dT/''t = -3£0T 4 /(cpR). (3) 

£ is the total hemispherical emissivity, a the Stefan-Boltzraann Constant, 
c the specific heat and p the density. In Fig. 1 the rate of cooling as a 
function of temperature is shown for a sphere of radius K = 0.5 cm with 
values of £ and pc typical of refractory metals. For the region of interest 
the rate varies quite widely, from 20 K/s to oearly iOUO ca it is 
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reasonable to consider measurements on samples aa snail as 1 on radius, the 
Inverse dependence on the radius, Eq. (3), would lead to cooling rates as 
hig h as 5000 K/s. The detectors will be operating at frequencies of about 
1 kHz, so the temperature will not be changing appreciably during the 
sampling period even for the highest cooling rate. Depending on sample 
size, the free-cooling period' will last between about 5 and 30 seconds for 
coolings to 2000 K and considerably longer for coolings extending below 
that temperature. 

The total power radiated by the sample, as shown In Fig. 1, gets quite 
large at the highest temperatures 'and will, of course, be eve., larger (or 
larger samples or ones with higher emissivities. At these levels, tin* 
heating of specimens to the desired temperature will become a major problem. 

For metals. Inductive heating will be the principal source, but to reach 
r.he highest temperatures it may have to be supplemented with a radiative 
source. Nonmetals will almost certainly require supplementary heating. 

The large amount of power being radiated by the sample, however, will 
simplify considerably the determination of that power at a distance. As a 
result, there should be no problem with signal level except possibly for 
the pyroelectric detectors with very small samples at the lowest temperatures. 

We foresee no insurmountable obstacle to the free-cooling experiments 
from the physical measurement aspect, that is, the radiation pyrometry 
measurement or the measurement of total radiant power with p roelectric 
detectors. Some of the problems in making these measurements have been 
pointed out in an earlier report ^ and are further discussed in this 
report. The real limitation to the free-cooling experiment* arise from the 
experimental approach. Because of the rapid cooling of the samples, temperature 
gradients will exist in the samples and for those of low thermal conductivity 


will be quite large. At some point these will mane it impossible to extract 
meaningful data from the measurements. We have spent considerable effort 
in the determination of the temperature gradients likely to be encountered 
in the experiments and what information can be extracted. The vaporization 
and/or decomposition of the sample trill determine the maximum temperatures 
which can be reached in these experiments. Both of these limitations on 
the applicability of the experiments will be discussed as well as the 
current plans for future development. 


Temperature Gradients 

The thermal gradients within a sphere, freely cooling by adiation, 
are given by the Fourier equation 


3 T 2 3T = 1 3T 
3i _2 r 3r a 3t 


(4) 


with the thermal diffusivity a = .'/pc, A being the thermal conductivity and 

pc the density times the specific heat, i.e., the heat capacity per unit 
volume. Solutions to Eq. (4) must satisfy the boundary condition 


^ 3 " ^4 

Ay r = eoT at r 


R, 


(5) 


where X 3T/3r is the heat conducted to the surface of the sphere which must 

equal the heat radiated from the surface given by the Stefan-Boltzmann 
A A 

relation eoT . The T term in the boundary condition, Eq. (5), makes the 

equation nonlinear and it cannot be solved analytically. Solutions to this 

[2 31 

problem have tc»en obtained by other workers ’ using the numerical method 
of finite differences, but their results are not available to us in a 
usable form. Ue are fortunate in having two generalized programs in the 
NBS computer library for the numerical solutions of nonlinear partial 
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differential equations. The programs are M0L1D (method of lines, one 
dimension) ^ and PDECOL (partial differential equations, collocation over 
piecewise polynominals) . The M0L1D program can be used to produce 
tabulated temperature values at various radii for a given set of times. Ue 
have used this program in generating the data discussed below. The program 
can be made very general in allowing for various initial temperature gradients 
and for different temperature dependences of the physical parameters, c, 
p, e, and X. The results are somewhat awkward to use in that one cannot 
conveniently generate data for different input parameters at the same 
temperatures. The PDECOL program produces a polynomial solution which may 
be better suited for this purpose but it has not been run to date. 

To describe and understand the temperature gradients in free-cooling 
spheres it has been extremely useful to define what we term the "limiting 
parabolic temperature distribution." Once a cooling sphere has this limiting 
distribution, its form will be maintained for the rest of the cooling cycle. 

To derive this limiting temperature distribution we assume first, that a gradient 
exists and second, that the gradients do not change as the 

sample cools. The second assumption is not strictly true, of course, and will be 

modified subsequently. In Fig. 2 we define the quantities describing the 

heat flow from a sphere; Q is the power radiating from the surface of 

the sphere of radius R, and q is the heat flux across a surface within 

the sphere at radius r. If all parts of the sphere are cooling at the same rate, 

and assuming that c, £, and X are all independent or only weak functions of 

T, then 


3 3 

4rrr _ r . 

q “ — 3 Q-“jQ* 
4tiR R 


( 6 ) 


This equation simply states that as the sphere cools, each volume element 
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gives up the sane amount of heat and that all the heat must be conducted to 
the surface at R. For q to be conducted across the surface at r, the 
temperature gradient at r must be such that 

q - -4irr 2 lg. . (7) 

Q Is given by the Stefau-Boltzman relation 

Q - 4 TrR 2 eoT g 4 ( 8 ) 

where T g indicates the surface temperature. 

Substituting (7) and (8) into (6) we have 

4 

eaT 

dT * - - jQ— rdr, ( 9 ) 

and integration from r to R yields 

eOT s 3 2 2 

T(r) - T,[l + ~2 Xr‘( R t >]- (10) 

The temperature at the center of the sphere relative to that at the surface is 


Usd 

T 

s 


EOT 3 R 


(id 


the temperature difference being 

4 

eaT R 

AT - T(o) - T g - — ~ ~ (12) 

These equations are accurate when the temperature difference, AT, is only a few 

percent of T , but not for larg.' values. Eq. (12) shows that the temperature 

difference is diminishing as T g drops whereas the initial assumption in Eq. (6) 

was that the difference was constant. To be correct, the temperature at the center 

of the sphere must drop more than the temperature at the surface. Therefore, 

q in Eq. (6) must be proportionately larger toward the center of the sphere 

by an amount which varies with r and T much as T(r)/T . Ve have fouad that 

8 

modifying Eq. (6) so that 
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9 r- 3 £a V ? ? ? r 3 

q - [T(r)/T s ] Z -*3 Q - [1 + (R -r ) J Q, (13) 

R R 

solving for a new T(r)/T using Eqs. (7) and ( 8 ), and by iterating 

8 

repeatedly, one obtains 

3 

T(r) - T g [l - — 2XR~< R ~ r ^) J ~ * (1A) 

which closely fits the coaputed results. 

In Fig. 3 we compare the temperature distributions derived from the partial 

differential equation (solid lines) with the predicted limiting parabolic 

distributions (dashed lines). «fe have assuned an initial condition of uniform 

temperature at 3000 K. The material parameters are X * 1.0 W/cm*K, 

3 

pc * 3.0 J/K*cm , and c * 0.3, values that correspond to tungsten and are 
relatively typical of most refractory metals. When the sphere initially 
starts to cool, the surface temperature drops rapidly until che temperature 
distribution within the sphere becomes parabolic and then the entire sphere 
commences to cool imiformly. This system reaches the limiting parabolic distribution 
in slightly over 0.5 seconds and after the temperature has fallen only about 2 1/2Z. 
It is interesting to note that the limiting temperature distribution depends only 
on the thermal conductivity and not on the heat capacity. Once the limiting 
distribution is established, the heat capacity determines the rate of coding only. 
The solution of the partial differential equation is needed to describe how the 
initial conditions die away and at what point the limiting parabolic behavior 
is applicable. 

Shows in Fig. 4 are results for spheres of three different radii for a 

hypothetical nometal where the material parameters, X » 0.046 W/cm*R, 

3 

pc - 4.6 J/K*cm , and £ ■ 0.3, are similar to those of solid alumina. The 
initial condition is again a itniform iperature of 3000 K and the material 
parameters are considered to be independent of temperature. In the figure 


94 



the time variable la reduced by dividing by che radius ot the spherical 

specimen (see Eq. (3)) which permits a bet :er comparison of the effect of 

the radius, rf c see that initially the suriai e temperature of the larger 

samples falls more rapidly (on our reduced time scale) until the limiting 

parabolic gradient is reached. The limiting distributions for the spheres are reached 

when t/R is approximately 3, 6, and 12 s/aa for the 0.125, 0.25, and 0.5 cm 

radii, respectively. At t/R ■ 20 the surface temperatures are nearly 

identical and the proportionality of the teaperature differences to the radii 

is obvious. That the three surface temperatures are the same at this time 

is only coincidental, for once the samples have achieved the llmltng 

parabolic distribution, the surface temperature of the smaller sphere is always 

cooling faster (in t/R units) than the larger sphere. This is a consequence 

of the initial rapid cooling of the surface of the larger sphere, so that 

it was losing proportionately less heat during subsequent periods than the 

mealier sphere. As a result, the average temperature of the larger sphere 

will always be greater than that of the smaller sphere. At lower temperaures, 

when the temperature gradients are small, the surface tasperature of the 

larger sphere will be higher than that of the smaller sphere. 

The magnitude of the change in the apparent heat capacities caused by 
the temperature gradients can be readily calculated when the gradients have 
the limiting parabolic form. The observed temperature change during a 
given time interval is that of the surface temperature, T^. The heat evolved 
during that interval, however, will depend on the change in the average 
temperature of the sample. The change in this average temperature, 

T ay , on cooling will always be greater than the change in T g because of the 
diminishing gradient with diminishing temperature. As a result, the apparent 
heat capacity, i.e., the observed energy radiated divided by the change in 
Tg, will always be greater than the actual heat capacity. The average 
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temperature for a sphere as determined by 


av 


1 2 f R 4 "r 2 T(r)dr » 


(15) 


(4/3)irR 


and using Eq. (14) for T(r) gives 

3T 


T 

av 


bR 


1 - 


tan 


-1 


— R 

a 


— R 
a 


(16) 


where a - 1 - eoT 3 R/(2A) and b = ea T 3 /(2AR). 161 
s s 

In Table 1, changes in T flv are given for each of two different radii, 
conductivities and surface temperatures. The results show that the increase 
in the apparent heat capacity over the actual varies, as might be expected. 


Table 1. Changes in T relative to changes in T p for spheres of two 

different radii and two different thernal conductivities including 
the percentage increase in the apparent heat capacity over the 
actual heat capacity 





R = 

0. 25 cm 


R * 

0.5 cm 


A(W/K*cm) 

T s (K) 

AT 

s 

T 

av 

AT 

av 

ZAC 

T 

av 

AT 

av 

ZAC 

0.1 

3000 

2950 

50 

3071. 24 
3016.50 

54.74 

9.5 

3147.58 

3087.49 

60.09 

20.2 

0.1 

2500 

2450 

50 

2533.87 

2481.20 

52.65 

5.3 

2569.09 

2513.57 

55.52 

n.o 

1.0 

3000 

2950 

50 

3006.91 

2956.46 

50.45 

0.9 

3013.87 

2962.96 

50.91 

1.8 

1.0 

2500 

2450 

50 

2503. 33 
2453.07 

50.26 

0.5 

2506.67 

2456.15 

50.52 

1.0 


linearly with R, inversely with A, and as the temperature cubed. The increases 
found for A - 0.1 W/K»cm are distressingly large for the conductivities of most 
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nonmetals will lie at or below this value. The conductivities of refractory 
aetals tend to lie between about 0.3 and 1.2 W/K*cm, so the Increase in 
apparent heat capacity will be sizable even for many metals. 

Because the heat capacity increase is nearly an exact linear function 
of the sample radius, it should be a relatively simple matter to correct 
for the increase by measuring samples of different size. Indeed, when the 
apparent heat capacity increase is large, it should be possible to determine 
thermal conductivities from a set of precise data. This process could be 
appreciably complicated by the temperature dependence of c, X, and 
e; these quantities were taken as constants in the above analysis, but they 
are usually relatively weak functions of T. This will work for solids but 
may break down completely for liquids if there is any stirring of the 
sample. If a sample is heated inductively, appreciable stirring action 
will occur during the heating period which may take time to die away. If a 
sample were radiantly heated in a low-gravity environment it may remain 
quiescent, but this will have to be determined. It is our intention, 
however, to use induction coils, operating at low power levels, to electro- 
magnetically position the samples within the apparatus. This action alone 
would probably add some stirring action to a liquid sample. 

Sample Vaporization 

In the high-temperature levitation calorimetry carried out by David 
Bonnell with John Margrave, ^ the upper temperature limit of the experiments 
vas determined by the vaporization of the refractory metal samples. In 
their case, the samples were in an inert-gas atmosphere and the vaporizing 
atoms condensing into a smoke which physically obscured the sample. This 
occurred at vapor pressures of about 10 Pa, which, in general, corresponds 
to temperatures about 200 K above the melting point of the metals. 
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The free-cooling experiments vill be in a vacuum so there will not be 
any smoke formation; sample vaporization, however, will still cause difficulties 
and determine the high temperature limit for the experiments. Assuming the 

evaporation coefficient for metals to be u * 1, the rate of evaporation at 

-3 2- 3 

10 Pa will be about W * 10 g/cm *s. Using a sample density of p = 12.5 g/cm , 

the rate of sample loss will be W/p * 8 x 10~^ cm/s. As a fraction of the 

total sample, this will be F = W4 ttR 2 / (4/3nR^p) = W/(Rp) = 8 x 10 ^/R s *. 

For R between 0.1 and 1 cm the loss will be less than O.lX/s and can probably 

be ignored. This vaporizing material will, however, be condensing on the 

apparatus. Consider a detector or detector objective lense 10 cm from a 

sample of R = 0.5 cm, the rate of deposition of the metal coating on these 

objects will be about 2 nm/s. This may seem a small amount but considering 

that an opaque metallic coating has a thickness of about 50 nm, indicates 

that the detectors will have to be protected if exnosed to samples vaporizing 

at the&: rates for periods of a few seconds or more. 

The heat loss from a sample will be rlgnificant at vapor pressures of 

10 Pa. Using AH v = 750 kJ/mol for the heat of vaporization of refractory 
[81 3 

metals, J and V =10 cm J /mol as the molar volume, the power of heat loss 

m 

2 

by vaporization at a vapor pressure o: 10 Pa is P v = WAlW (pV^) * 6 W/cm . 

In Table 2, P y is compared with the radiative power, P R , for 3urfj»''«-.s at 
three different temperatures. 

Table 2. Power of evaporative heat loss (vapor pressure * 10 Pa) compared 
with radiative power (i • 0.3) at various sutface temperatures 


T(K) 

P (W/cm 2 ) 
v 

Pr(W/ cm 2 ) 

P y /P R x 100 (Z) 

2000 

27 

6 

22 

3000 

137 

6 

4.4 

4000 

435 

6 

1.4 


i 
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It is apparent that with samples at temperatures where the vapor pressure 
is 10 Pa, a significant correction will have to be made for sample . 
zation, and for those sample materials reaching that pressure at ‘ w 
temperatures, the correction will have to be of fairly high accuracy. 

There are also potential problems from the sample v.»por due to its 
absorption spectrum. Nearly all materials will have absorption bands in 
the region of interest, 0.3 pm to 15 pm for the total radiance measurement . 

The atomic absorption lines of metal vapors will probably remain sharp up to 
pressure of 10 Pa and block out very little of the total radiation. Above 
that pressure, however, '■.he lines will begin to broaden and an appreciable 
amount of radiation from the sample would be blocked from the detector. For 
the radiation pyrometry, which will use narrow bandwidths between 0.4 and 
1 pm, care will have to be taken that absorption lines do not occur at the 
wavelengths being used. 

This discussion on vaporization has been restricted to metals because of 
the availability of data and the generalizations that can be made. The same 
type of problems and limitations will occur with nonmet -.111c samples and will 
probably be much more involved and stringent. The evaporation coefficient will 
most likely be less than 1 and vaporization can be accompanied by dissociation 
of the molecular species so that quantities such as rates and heats of 
vaporization are less well defined. Molecular species in the vapor will make 
the absorption spectrum much more complex with numerous broad absorption bands. 
Finally, if decomposition is occurring at the surface, the radiative properties 
of the surface could be altered considerably and introduce errors in interpre- 
tation of the radiation measurements. 
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Total Radiance Measurement 


The unique aspect of the free-cooling experiment and that with least 
precedent, is the proposal to use pyroelectric detectors to measure the 
total radiant heat loss from the sample. The detectors must be uniformly 
black over the wavelength 0.4 to 14 im which will include 99% of the radiation 
from samples at tanperatures between 1500 and 4000 K. When the field of 
view of the detector encompasses the entire sphere of the sample, the 
signal will be proportional to the total hemispherical emittance, ^ the 
proportionality will be the area of the detector relative to the area of 
the sphere whose radius is the distance from the detector to the center of 
the sample. 

A pyroelectric is a poled ferroelectric, which, when heated, changes 
its spontaneous polarization. LiTaO^ detectors \ 11 be used because the 
883 K Curie temperature of this material makes it well suited for a hi^i, 
temperature apparatus and they are available in a wide range of produces 
from several manufacturers. A pyroelectric detector has electrodes on two 
sides of a thin wafer of material forming a capacitor, a change in polarization 
of the material can then be used to generate a current in an external 
circuit. The current is proportional to the rate of change of temperature 
so that radiation from a continuous source must be chopped to get useable 
response. A black coating is applied to the front surface of the pyroelectric 
to enhance the light absorption and the back surface is bonded to the 
mounting which serves as a heat sink. Fig. 5 shows the thermal response of 
an idealized pyroelectric detector. The diagram at the top of the figure 
describes the model used; P is the power of the incident radiation, c is 
the thermal capacity of the pyroelectric which is at temperature T, h is 
the thermal conductance to the heat sink at T . The response of the model 
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is given by 


P - h(T-T ) 
o 

P - h(T i -T o ) 


e c 


( 1 ?) 


which is plotted as curve A in Fig. 5 with P 0 and the initial temperature, 

T, > T and a3 curve B with P * P and T. ** T . The thermal response of the detector 
i o i o r 

to chopped radiation is represented ly the saw-toothed curve which rises from T when 
the radiation is turned on until reaching a steady state. At the steady 
state, the temperature of the pyroelectric will be following the path 
indicated by the arrows at the intersection of curves A and R in Fig. 5. 

For this idealized example, che electrical response will be a square wave 
when the chopping period is sufficiently small compared with the thermal 
relaxation of the device, T = c/h. As the chopping period is made longer, 
the curvature in the thermal response will became more evident and the 
square wave will begin to droop. 

The model for the thermal response is idealized in several ways. 

Foremost is that the model implies that the pyroelectric is uniformly 
heated whereas it actually responds to the thermal pulse as it propagates 
through the material. Also, the black coating introduces a thermal lag 
whicn together with electrical response of the circuit produce a sizable 
rise-time for the device. Fig. 6 shows oscilloscope tracings of the response 
of a pyroelectric detector at four different chopping frequencies. At 300 
Hz the rise time of the response is very evident but part of this is due to 
the transition period of the chopper. The constant voltage period is 
somewhat less than a millisecond. At 150 Hz the voltage is essentially 
constant for more than 11/? m3. At 75 Hz the droop it che voltage signal 
is evident and from the 15 Hz curve the thermal time constant for the 
device of 15-20 ms can be dedu- .d. It was intended that this particular 
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detector could be used et 300 Hz with the voltage sampled for 1 as during 
each half 'cycle, the 300 Hz curve in Fig. 6 shows that either the sampling 
period or the chopping frequency will have to be reduced to get the full 
signal. 

Gold black is used for the absorptive coating on the pyroelectric 
detectors and numerous reflectance measurements have shown these coatings 
to be uniformly black to between 98 and 100Z for wavelengths from 0.3 to 
10 pm. Beginning between 10 and 20 pm there is a drop-off to values as low 
as 90Z with the position of the drop-off varying from sample to sample. At 
wavelengths above 10 p, however, we are concerned only with the long wavelength 
tail at the lowest temperature so the errors introduced should be small. 

Because of the smallness of the detectors and the degree of difficulty in 
making accurate reflectance measurements, the blackness of individual 
detectors cannot be determined directly and must be assumed to be similar 
to the measured coatings. 

The calibration of the pyroelectric detectors will present some diffi- 
culties. The response of the detectors are supposed to be linear over the 
range of power to be covered by the experiments for samples of 0.1 to 1 cm 
radius and temperatures from 1500 to 4000 K. The linearity and magnitude 
can be established to a fair degree of precision using standard tungsten 
filament lamps. There will be, however, uncertainties of several percent 
arising from the transfer of the calibration as well as the original lamp 
calibration itself. 

Current Experimental Pl an 

Because of the untried nature of the pyroelectric detection system in 
particular and the free-cooling experiment in general, we are proceeding 
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with a ground baaed simulation of the apace experiment. The setup to be 
used is shown in Fig. 7 where in place of the sample being levitated within 
the apparatus, it is simply suspended from a wire of the sme material. 

Tungsten is the material of choice since with its high melting point, 

3650 K., experiments could probably be run at temperatures up to 3500 K or, 
as is more likely, up to the limit of the available power for Induction 
heating. The heat capacity and emittance are reasonably well known for 
solid tungsten so that the observed results for the free-cooling of the 
saaple can be compared with calculated results. The correction for the 
heat leak to the suspension wire will probably be relatively small and, in 
any case, it can be minimized by electrically heating the wire. A hole 
drilled in the sample can .e used for blaclcbody temperature measurements so 
that it will not be necessary to rely on any assimptions about the surface 
emittance in order to determine the sample temperature. The simulated 
experiments will allow investigation of many of the problem areas such as 
reflected light getting into the wide aperture of the pyroelectric detector 
and light reflected from the coils back onto the sample. 

A high-speed, 12 bit, digital voltmeter will be used to record the amplified 
pyroelectric detector signals. The unit is controlled and the data logged 
by a desfr-top computer. Pyrametry will initially be carried out with an 
existing millisecond pyrometer, but we are about to start designing a new 
pyrometer which will be dedicated to the free-cooling experiments. The 
pyrometer will be a two or three color instrument sampling at millisecond 
Intervals. Control and data logg ing will be handled by the same system 
used for the pyroelectric detector. 
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Conclusions and Discussion 


The upper temperature limit of the free-cooling experiments trill be 
limited by the vaporization of the sample. This trill occur when the vapor 
pressure of the material reaches the vicinity of 10 Pa. 

For nonmetals there will be problems from numerous sources. The low 
electrical conductivity may preclude inductive heating and limit the ability 
to use radio frequency fields for positioning the samples. The transparency 
of some nonmetallic materials trill cause difficulties in temperature 
determinations^. Temperature gradients, determined by the thermal conduc- 
tivity of the material, will limit the precision of the derived physical 
quantities. For A < 0.01 V/K*cm, it ray be difficult to extract any 
meaningful results from the data. For 0.01 < A < 0.1 W/K*cm, it should be 
possible to obtain both heat capacities and thermal conductivities of 
moderate precision for solid samples. For X > 0.1 W/K*cm, relatively 
precise heat capacities should be obtained together with a low precision 
value of A, again for solid samples. Meaningful data on liquid samples 
with X much below 1 W/K*cm could only be obtained in the limits of rapid 
stirring or no stirring. It will probably not be possible to derive heats 
of fusion from the cooling curves of poorly conducting materials. 

Nith metals, vaporization of the sample appears to be the only impediment 
to heat capacity measurements by the free-cooling technique. The high 
thermal conductivity of the samples will keep the gradients small enough 
that stirring within liquid samples will produce only a minor perturbation 
on the measurements. It is unfortunate that with few exceptions (Hf, Zr, 
and Ft), the vapor pressure of high-melting metals will limit the experiments 
to temperatures just slightly above their melting points, because it would 
be of great vj-lue to obtain the heat capacities of liquid metals over a 
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wide temperature range. The free-cooling experiments do, however, open a 
new prospect; to measure the heat capacities of the liquid metals into the 

f Q] 

supercooled range. Container less samples should, and have been observed 1 
to, supercool to a considerable extent. The free-cooling experiment will 
permit the heat capacity to be determined all the way to' the point where 
solidification commences and the sample recalesces. It will also be possible 
to change the sample size and thus the rate of cooling, which may produce 
variations in the extent of supercooling and even possibly the heat capacity. 
The prospect of the supercooling experiments with metals makes the development 
of the free-cooling technique worthwhile even if the determination of 
nonmetal heat capacities were to turn out to be infeasible. 
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Fig. 3. Temperature gradients in a sphere of 1 ui radius freely cooling by 
radiation from an initial constant temperature of 3000 K. The 
material parameters, e - 0.3, pc - 3.0 J/K'cm 3 , and X - lW/K*cm 
correspond closely to those of tuigsten. The solid curves are 
solutions of the partial differential equation describing the system. 
The dashed curves are the limiting parabolid temperature 
distributions for the ame surface emperat ures. 
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Temperature gradients in freely cooling spheres of three different 
radii of a hypothetical material with e ■ 0.3, pc - 4.6 J/K*cia 3 and 
X “ 0.046 rf/K*<m. The initial condition is a uni fore temperature 
of 3000 K. The dash-dot line at the bottom is a reference line to 
highlight the direct proportionality of the magnitude of the 
gradients to the radii. 








Thermal behavior of an idaalisad pyroelectric detector 





Pig. 6. Oscillograph tracing of the output of a pyroelectric detector at four 
different light chopping frequencies. 
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Fig. 7. Diagram of simulated space free-cooling experiment. 
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